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Quantum mechanics and the theory of gravity are presently not compatible. A particular question 
is whether gravity causes decoherence - an unavoidable source of noise. Several models for gravi¬ 
tational decoherence have been proposed, not all of which can be described quantum mechanically 
[1-3]. In parallel, several experiments have been proposed to test some of these models [4-6], where 
the data obtained by such experiments is analyzed assuming quantum mechanics. Since we may 
need to modify quantum mechanics to account for gravity, however, one may question the validity of 
using quantum mechanics as a calculational tool to draw conclusions from experiments concerning 
gravity. 

Here we propose an experiment to estimate gravitational decoherence whose conclusions hold even 
if quantum mechanics would need to be modified. We first establish a general information-theoretic 
notion of decoherence which reduces to the standard measure within quantum mechanics. Second, 
drawing on ideas from quantum information, we propose a very general experiment that allows 
us to obtain a quantitative estimate of decoherence of any physical process for any physical theory 
satisfying only very mild conditions. Finally, we propose a concrete experiment using optomechanics 
to estimate gravitational decoherence in any such theory, including quantum mechanics as a special 
case. 

Our work raises the interesting question whether other properties of nature could similarly be 
established from experimental observations alone - that is, without already having a rather well 
formed theory of nature like quantum mechanics to make sense of experimental data. 


Experiments aiming to test the presence - and amount 
- of gravitational decoherence generally go beyond estab¬ 
lished theory. Many theoretical models for gravitational 
decoherence have been proposed [1, 3, 7-15], and it is 
wide open if one of these proposals is correct. As such, 
experiments are of a highly exploratory nature, aiming 
to establish data points to which one may tailor future 
theoretical proposals. This task is made even more dif¬ 
ficult by the fact that quantum mechanics and gravity 
do not go hand in hand, and indeed quantum mechan¬ 
ics may need to be modified in a yet unknown way in 
order to account for gravitational effects such as deco¬ 
herence. We are thus compelled to design an experiment 
that provides a guiding light for the search for the right 
theoretical model - or indeed new physical theory - whose 
conclusions do not rely on quantum mechanics. 


Decoherence in QM 

As an easy warmup, let us first focus on the concept of 
decoherence within quantum mechanics. We first show 
how the protocol given in Figure 1 allows us to esti¬ 
mate quantum mechanical decoherence without knowing 
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the decoherence process, and without doing quantum to¬ 
mography to determine it. Traditionally, the presence 
of decoherence within quantum mechanics is related to 
the change of state due to measurement and the ’’col¬ 
lapse of the wavefunction”. There are two complimentary 
ways to view this based on unconditional and conditional 
states. Given some pure quantum state, cr|0) +/?11), and 
an arbitrarily accurate measurement of the variable di¬ 
agonal in this basis, the post-measurement conditional 
states are |0) or |1) conditioned on the measurement 
outcome. On the other hand if this measurement has 
taken place but the results are unknown, the resulting 
unconditional state is given by the a quantum density 
matrix p = |a| 2 |0)(0| + |/3| 2 |1)(1|. The vanishing of the 
off-diagonal matrix elements in the measurement basis 
for the post measurement unconditional state forms de¬ 
coherence. If the measurement is not arbitrarily accurate 
(i.e. weak) the off-diagonal matrix elements are reduced 
but do not vanish. More general forms of decoherence 
correspond to a decay of off-diagonal terms in the den¬ 
sity operator p with respect to any basis, and can occur 
due to the interaction of the system with an environment 
that may not be a measurement procedure of any kind. 
It is clear that this way of thinking about decoherence is 
entirely tied to the quantum mechanical matrix formal¬ 
ism, and also offers little in the way of quantifying the 
amount of decoherence in an operationally meaningful 
way. 
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FIG. 1: Any decoherence process - also known as a (quan¬ 
tum) channel - can be thought of as an interaction Ui of 
the system A' with an environment Ei n . In quantum me¬ 
chanics, the resulting state is the output of the channel 
p B = r a'^b(Pa') = Tr E [Uip AI ® \U\]. In gen¬ 

eral, B may be a smaller or larger system than A' . In the 
examples below, however, we will focus on the case where 
A' and B 1 have the same dimension, corresponding to the 
case where a fixed system A' = B experiences some inter¬ 
action with another system E ln = E. The channels abil¬ 
ity to preserve quantum information - that is, the amount 
of decoherence - can be characterized by how well it pre¬ 
serves entanglement between an outside system A and A'. 
We note that our treatment of theories that go beyond stan¬ 
dard quantum mechanics makes no statement whether the 
environment is an actual physical system, or merely a mathe¬ 
matical Gedankenexperiment possibly used to describe an in¬ 
trinsic decoherence process. In full generality, the experiment 
consists of a Bell experiment in which a source of decoher¬ 
ence is introduced deliberately. For simplictly, we consider 
an experiment for the CHSH inequality, although our analy¬ 
sis could easily be extended to any other Bell inequality. In 
each run, a source prepares the maximally entangled state 
&aa', where A' is subsequently exposed to the decoherence 
process to be tested. We then perform the standard CHSH 
measurements: system A is measured with probability 1/2 
using observables Ao = ax and Ai = az respectively. Sys¬ 
tem B is measured using observables Bo = (ax — az)l\/^ 
and B\ = (ax + az)/\/2 with probability 1/2 each. Per¬ 
forming the experiment many times allows an estimate of 
13 = Tr [pab(A 0 ® Bo + A 0 0 Hi + Ai 0 Bq - Ai 0 Bi)]. 


The modern way of understanding decoherence in 
quantum mechanics in a quantitative way is provided 
by quantum information theory. One thereby thinks of 
a decoherence process as a interaction of a system A' 
with an environment as described in Figure 1, resulting 
in a quantum channel Ta'^b- The amount of of deco¬ 
herence can now be quantified by the channel’s ability 
to transmit quantum information, i.e., its quantum ca¬ 
pacity (see e.g. [16] and the appendix for further back¬ 
ground) . Concretely, one considers n uses of the channel 
given by and asks how many qubits k we can 

send in relation to n using an error-correcting encoding. 
Of particular interest is thereby the so-called single shot 
capacity, which determines the largest rate k/n up to a 


given error parameter for any choice of n b This single¬ 
shot capacity is determined by the so-called min-entropy 
H min (A|£;) [17, 18J. 

Apart from its information-theoretic significance, the 
min-entropy has a beautiful operational interpretation 
that also makes its role as a decoherence measure in¬ 
tuitively apparent. Very roughly, the amount of deco¬ 
herence can be understood as a measure of how cor¬ 
related E becomes with A. Suppose we start with a 
maximally entangled test state &aa' where the decoher¬ 
ence process is applied to A!. This results in a state 
1*1 ’ abe) (see Figure 1). If no decoherence occurs, the 
output state will be of the form <&ab |0)(0|e where 
A’ = B. That is, A and B are maximally entangled, but 
A and E are completely uncorrelated. The strongest de¬ 
coherence, however, produces an output state of the form 
*&AE\ ® Pe 2 <S> |0)(0|b where A! = E\ and E = EiE 2 . 
That is, A is now maximally entangled with E\, whereas 
A and B are completely uncorrelated. What about the 
intermediary regime? The min-entropy can be written as 

H min (A|£) = - log d A Dec(A|F;) , (1) 

where dA is the dimension of A, and [19] 

Dec(A\E) = max F 2 ($ A aa It ® U e ^a'(pae)) ■ (2) 

The maximization above is taken over all quantum oper¬ 
ations IZe^A' on the system E, which aim to bring the 
state pae as close as possible to the maximally entan¬ 
gled state <5>aa'- Intuitively, Dec(A|i5) can thus be un¬ 
derstood as a measure of how far the output pae is from 
the setting of maximum decoherence (where pae = $ ae 
is the maximally entangled state). If there is no decoher¬ 
ence, we have pae = l/dA®p_E giving Dec(A|£’) = 1 /d 2 A 
and H m ; n (A|£ l ) = logd^- If there is maximum decoher¬ 
ence, we have pae x = <&AA' giving Dec(A|£’) = 1 and 
H m i n (A|i?) = — log dA where TI e ^a' = Tr B2 is simply 
the operation that discards the remainder of the environ¬ 
ment E 2 . A larger value of Dec(A|i?) thus corresponds 
to a larger amount of decoherence. In the quantum case, 
T3ec(A\E) can be computed using any semi-definite pro¬ 
gramming solver [20, 21]. 

We hence see that in quantum mechanics, the rele¬ 
vant measure of decoherence is simply Dec(A|.E). How 
can we estimate it an experiment? Our goal in deriv¬ 
ing this estimate will be to rely on concepts that we can 
later extend beyond the realm of quantum theory, deriv¬ 
ing a universally valid test. It is clear that to estimate 
Dec(A|.E) we need to make a statement about the entan¬ 
glement between A and E - yet E is inaccesible to our 
experiment. A property of quantum mechanics known as 
the monogamy of entanglement [22] nevertheless allows 


1 The asymptotic quantity often considered in information theory 

arises in the limit of n —> oo, but the single-shot capacity gives 
more refined statements which are valid for any n. 
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such an estimate: if pab is highly entangled, then pae 
is necessarily far from highly entangled. Since low en¬ 
tanglement in pae means that Dec(A|.E) is low, a test 
that is able to detect entanglement between A and B 
should help us bound Dec(A\E) from above. We note 
that whereas any experimental proposal demands that 
we specifiy concrete measurements to be performed, our 
conclusions remain valid even if we do not have full con¬ 
trol over the measurements, possibly because they are 
also somehow affected by an gravitational interaction in 
an unknown fashion. Dealing with unknown states and 
measurements is the essence of so-called device indepen¬ 
dence [23] in quantum cryptography. Allowing arbitrary 
measurements again forms a crucial stepping stone, en¬ 
abling us to extend our results beyond quantum mechan¬ 
ics. 


Beyond QM 

The real challenge is to show that the conclusions of 
our test remain valid even outside of quantum mechanics. 
Since we want to make as few assumptions as possible, we 
consider the most general probabilistic theory, in which 
we are only given a set of possible states 12 and measure¬ 
ments on these states. Every measurement is thereby a 
collection of effects e a : 12 —>• [0,1] satisfying e a (u) ^ 0 
and J2 a e a (oj) = 1 for all w £ 12. We also refer to a 
measurement as an instrument Ad = {e a } a - The label a 
corresponds to a measurement outcome ’a’. The notion 
of separated systems A, B and E is in general difficult to 
define uniquely. We thus again make the most minimal 
assumption possible in which we identify ’’systems” A, B 
and E by sets of measurements that can be performed. 
For simplicitly, we take measurements and operations in 
the sets A,B, and E to commute, but do not impose any 
other strucuture. We thus merely use labels A and B 
and E for commuting measurements. This means that 
for maps going from a system E to an output system A' 
like TZe^A' the map is really from E to E and we use A' 
merely to remind ourselves we consider a restricted class 
of measurements on the output. Again, this is analogous 
to quantum mechanics where such measurements consist 
of operators on A' and the identity elsewhere (see ap¬ 
pendix for a discussion). 

The first obstacle consists of defining a general no¬ 
tion of decoherence. We saw that quantumly decoher¬ 
ence can be quantified by how well correlations between 
A and A' are preserved, and this can be measured by 
how well the decoherence process preserves the maxi¬ 
mally correlated state. Indeed, we can also quantify clas¬ 
sical noise in terms of how well it preserves correlations, 
where the maximally correlated state takes on the form 
(1 /gU) J2 a l a )(°U ® | a)(a|A' for some classical symbols 
a. We hence start by defining the set of maximally cor¬ 
related states, by observing a crucial and indeed defining 
property of the maximally correlated in quantum me¬ 
chanics. Concretely, A and A' are maximally entangled 


if and only if for any von Neumann measurement on A, 
there exists a corresponding measurement on A' giving 
the same outcome. Again, the same is also true classically 
but made trivial by the fact that only one measurement is 
allowed. In analogy, we thus define the set of maximally 
correlated states as 

i>AA' = G n AA ’ | VAT 4 = (e 4 } a 3M b = {ef } 0 

such that ^e^ef(w) = 1 j . (3) 

This set coincides with the set of maximally entangled 
states in quantum mechanics, where A! can potentially 
contain an additional component a a' in <&aa'®&A' which 
is irrelevant to our discussion. We thus define 

Dec(A|£% = sup sup F 2 ($ A A', He^A'{uae)) , 

( 4 ) 

where loae is the state shared between A and E according 
to the general physical theory. The fidelity between two 
states wi and u >2 is thereby defined in full analogy to the 
quantum case [25] as 

F(wi,W 2 ) = inf , (5) 

M 

where the minimization is taken over all possible mea¬ 
surements M, and M(u) denotes the probability distri¬ 
bution over the measurement outcomes of Ad. That is, 
the fidelity can be expressed as the minimum fidelity be¬ 
tween probability distributions of classical measurement 
outcomes. How about the transformation IZe^A’ 2 In 
general, it is difficult to characterize the set of allowed 
transformations 7 Ze^a' in arbitrary physical theories, 
however we will not need make IZe^a' explicit in or¬ 
der to bound Dec(A\E). Equation (4) gives us the fa¬ 
miliar quantity within quantum mechanics, but provides 
us with an a very intuitive way to quantify decoherence 
in any physical theory that admits maximally correlated 
states. We emphasize that with our general techniques 
the latter demand could be weakend to allow all theories, 
even those who only have (weak) approximations of max¬ 
imally correlated states. However, as we are not aware 
of any physically motivated example of such a theory, 
we leave such an extension to future study for clarity of 
exposition. 

The second challenge is to prove that our test actu¬ 
ally provides a bound on Dec(A|E) w . Note that without 
quantum mechanics to guide us, all that we could rea¬ 
sonably establish by performing measurements on A and 
B are the probabilities of outcomes a and b given mea¬ 
surement settings x and y. That is, the probability 

Pr[a, b\x, y ] w = efef ( uj A b ) , (6) 

where G M£ and ef G M B . Yet, given the sys¬ 
tem E is entirely inaccessible to us we have no hope of 
measuring Pr[o, b, c\x, y, z] u directly, where z denotes a 
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FIG. 2: Comparison of the quantum bound with the actual values of Dec(A\E) p for some example channels 
and measurements (colors online). The black dashdotted line on top shows the quantum bound, i.e. the maximal value of 
Dec(d|f?) p that is compatible with a measured CHSH value in quantum theory. The other four plots are parametric plots: 
The parameter that is varied is the noise parameter of the channel. For each noise parameter, the value of Dec(d|f?) p of the 
resulting state is calculated, as well as the CHSH value /3 that one would meausre for this state using the standard measurements 
in the X-Z -plane that would be optimal for an EPR pair. This measurement happens to be optimal for the resulting state for 
the depolarizing channel, but not for the dephasing channels. The orange solid line also shows such a parametric plot for the 
dephasing channel, but for that line, the CHSH value /3 is not calculated for the standard measurement for the EPR pair but 
for the measurement that is optimal for the actual resulting state. This is done using a formula found in [24]. The resulting 
curve is independent of the dephasing direction. 


measurement setting on E with outcome c. Neverthe¬ 
less, similar to quantum entanglement, it is known that 
no-signalling distributions are again monogamous [26] - 
and it is this fact that allows us to draw conclusions about 
E by measuring only A and B. We will therefore make a 
non-trivial assumption about the physical theory, namely 
that no-signalling holds between A, B and E. We em¬ 
phasize weaker constraints on the amount of signalling 
could also lead to a bound - but we are not aware of any 
other concrete example to consider. Mathematically, no 
signalling means that the marginal distributions obey 

Va, x, y, y', z, z' Pr[a|x, y, z ] w = Pr[a|x, y', z'] u , (7) 

that is, the choice of measurement settings y, y' and z, z' 
does not influence the distribution of the outcomes a. A 
set of distributions is no-signalling if such conditions hold 
for all marginal distributions. 

Abstract experiment 

Our method is fully general and can in principle be 
used to measure the decoherence of any physical process. 
Figure 1 illustrates the general procedure. We create 
an entangled pair, and use half of this entangled pair 
to probe the unknown decoherence process. To estimate 
Dec{A\E) we will make use of the fact that in QM entan¬ 


glement is monogamous, or more generally - when con¬ 
sidering theories beyond QM - that no-signalling correla¬ 
tions are monogamous. This allow us to make statements 
about the correlations between A and E, even though we 
can only perform measurements on A and B. A test that 
allows us to bound Dec(A|.E) from observations made 
on A and B alone is given by a Bell inequality [31, 32]. 
For the purpose of illustration, we consider creating an 
entangled state <&aa' and perform a test based on the 
CHSH inequaltity [33] (see Figure 1). We emphasize that 
our methods are fully general and could be used in con¬ 
junction with other inequalities and higher dimensional 
entangled states. 

As an easy warmup, let us first again consider what 
happens in quantum mechanics. For now, we assume 
that the measurement devices have no memory. That is, 
the experiment behaves the same in each round, inde¬ 
pendent on the previous measurements. It is relatively 
straight forward to obtain an upper bound on Dec(A|A') 
by extending techniques from quantum key distribution 
(QKD) [23]. In essence, we maximize Dec(A|A') over all 
states that are consistent with the observed CHSH cor¬ 
relator fi (see Figure 1). This maximization problem is 
simplified by the inherent symmetries of the CHSH in¬ 
equality, allowing us to reduce this optimization problem 
to consider only states that are diagonal in the Bell ba¬ 
sis. We proceed to establish properties of min and max 
entropies for Bell diagonal states, leading to an upper 
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(a) QM region, GPT region and forbidden region 



(b) Plot of the GPT bound on Dec(T|I?)aj 


FIG. 3: Allowed values of the decoherence quantity for 
measured CHSH values (colors online). Part (a) shows 
what values of the decoherence quantity are compatible with 
some measured CHSH value /3, assuming either quantum the¬ 
ory or any other probabilistic theory. The dark green region 
consists of all points (f3, Dec(A|_E) p ) for which there exists 
a quantum state pab and two pairs (Ao,Ai) and ( Bq,B \) 
of observables with the according values, i.e. the bound is 
tight. The red region shows pairs (/3,Dec(A\E) UJ ) that can¬ 
not be realized in any non-signalling probabilistic theory, but 
not necessarily all of them. In other words, the curve between 
the light green area and the red area is a bound on Dec(A|l?)„ 
which is valid for all non-signalling probabilistic theories, but 
the light green area might still contain forbidden pairs. What 
is certain is that a bound for any specific non-signalling prob¬ 
abilistic theory runs below the red region. Part (b) shows 
a zoomed-in plot of the border line between the forbidden 
region and the region which is potentially allowed by GPTs, 
making it easier to see where the line runs for values of /3 close 
to 2. In a world constrained only by no-signalling, 8 = 4 is 
possible [27-30], 


bound. Concretely, we show in the Appendix (Theo¬ 
rem B.l) that 


Bec(A\E) < h{P) , (8) 

where h(/3) is an easy optimization problem that can be 
solved using Lagrange multipliers. We have chosen not to 


weaken this bound by an analytical bound that is strictly 
larger, as it is indeed easily evaluated (see Figure 3). If 
the devices are allowed memory, then a variant of this 
test and some more sophisticated techniques from QKD 
nevertheless can nevertheless be shown to give a bound. 



FIG. 4: Probing an optomechanical system. Our goal 
is to create two entanglement between two opto-mechanical 
cavities. One cavitiy thereby has a movable mirror that intro¬ 
duces gravitational decoherence. Two cavities each contain a 
Raman single photon source controlled by an external laser 
‘write field’ E(t) [34], This write-field is used to map excita¬ 
tions in the atomic sources to single photon excitations in the 
cavities. The top cavity has fixed end mirrors while the bot¬ 
tom cavity has one mirror that is harmonically bound along 
the cavity axis and can move in response to the radiation 
pressure force of light in the cavity. The Raman sources are 
first prepared in an entangled state. This setup is a modifi¬ 
cation of the one proposed by Bouwmeester [5] in which an 
itinerant single photon pulse is injected into a cavity rather 
than created intra-cavity as here. Our modification avoids the 
problem that the time over which the photons interact with 
the mechanical element is stochastic and determined by the 
random times at which the photons enter and exit the cavity 
through an end mirror. In the new scheme, the cavities are 
assumed to have almost perfect mirrors — very narrow line 
width [35] (see appendix for details). 


How can we hope to attain an estimate outside of quan¬ 
tum mechanics? Let us first give a very loose intuition, 
why performing a Bell experiment on A and B , may al¬ 
lows us to bound Dec(A|£i) w . It is well known [26] that 
non-signalling correlations are also monogamous. That 
is, if we observe a violation of the CHSH inequality as 
captured by the measured parameter /3, then we know 
that the violation between A and E and also between E 
and B must be low. Note that the expectation values 
Tr [pab(A x ® By)] in terms of quantum observables A x 
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and By can be expresssed in terms of probabilities as 
Tr [pab{A x <S> By)} 

= ^2 Pl [ a , a \ x ,y]u ~ Pr[a, ~a\x,y] u , (9) 

a£{±l} 

where we have again used ujab in place of pab to remind 
ourselves that we may be outside of QM. In fact, if (3 is 
larger than what a classical theory allows (/? > 2), then E 
and B cannot violate the CHSH inequality at all. Let us 
now assume by contradiction that the state ujae shared 
between A and E would be close to maximally correlated. 
Then by definition of the maximally correlated state, for 
every measurement on A, there exists some measurement 
on E which yields (almost) the same outcome. Hence, if 
ujae would be close to maximally correlated, then we 
would expect that E and B can achieve a similar CHSH 
violation than A and B - because E can make measure¬ 
ments that reproduce the same correlations that A can 
achieve with B. Yet, we know that this cannot be since 
CHSH correlations are monogamous. 

In the appendix, we make this rough intuition precise. 
While we do not follow the steps suggested by this intu¬ 
ition, we employ a technique that has also been used for 
studying monogamy of CHSH correlations [26]. Specif¬ 
ically, we use linear programming as a technique to ob¬ 
tain bounds. We thereby first relate the fidelity to the 
statistical distance, which is a linear functional. We are 
then able to optimize this linear functional over proba¬ 
bility distributions Pr[a, b, c\x, y, z] u satisfying linear con¬ 
straints. The first such constraint is given by the fact that 
we consider only no-signalling distribtions. The second 
by the fact that the marginal distribution Pr[a, b\x, y]ai 
leads to the observed Bell violation /3. The last one stems 
from the fact that maximal correlations can also be ex¬ 
pressed using a linear constraint. Solving this linear pro¬ 
gram for an observed violation /3 leads to Figure 3. 


Optomechanics experiment 

To gain insights into the significance of gravitational 
decoherence, we examine Diosi’s theory of gravitational 


decoherence [1]. This is equivalent to the decoherence 
model introduced in Kafri et al. [9]. It can be applied 
to an optomechanical cavity in which one mirror is free 
to move in a harmonic potential with frequency uj m as 
in Figure 4. The master equation for a massive particle 
moving in a harmonic potential, including gravitational 
decoherence is 

^ = -iu> m [b'b, p\ - A [b + b\[b + b\p\] (10) 


where 


mujm 


1 


h = \l + i . -p (11) 

V Zn v Zrimujm 

with x,p the usual canonical position and momentum 
operators for the moving mirror. We have that 


A = A„ 


- A 


heat ? 


( 12 ) 


where the gravitational decoherence rate A grav is given 

by 


A 


grav 


2t r G A 
3 uj m 


(13) 


with G the Newton gravitational constant and A the den¬ 
sity of the moving mirror. As one might expect A grav is 
quite small, of the order of 10~ 8 s” 1 for suspended mir¬ 
rors with u] m ~ 1. The term 


Aheat hQ , (14) 

with Q = uj/^m corresponds to mechanical heating. To 
see effect of the gravitational term stand out next to the 
mechanical heating we thus need to make the tempera¬ 
ture T low. A calculation shows that this model leads to 
a dephasing channel T(p) = pp + (1 — p)ZpZ 1 where p is 
a function of the density A, and the other parameters. 
In the appendix, we show that for this model 


Dec( A\E) P = - 1 + * 1 — exp (-4(1 + 2 



(15) 


where G is the Newton gravitational constant, ks is the 
Boltzman constant, and h the Planck constant (see Fig¬ 
ure 5 for the other parameters), (see Figure 5 for param¬ 
eters) 


Discussion 

What have we actually learned when performing such 
an experiment? We first observe that the measured /3 
always gives an upper bound on the amount of decoher- 
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T = 1 nK 
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T = 5 nK 


T = 10 nK 


T = 50 nK 





t t t 


- Gravitational decoherence included, density A = 2,102 X 10 4 kgm 3 (rhenium) 

-Gravitational decoherence included, density A = 2,7 x 10 3 kgm~ 3 (aluminum) 

-Gravitational decoherence neglected, (material-independent mechanical heating only) 


FIG. 5: Predicted values of the decoherence quantity in the optomechanical experiment. This figure shows the 
predicted values of Dec(A\E) p as a function of the running time of the optomechanical experiment for different temperatures 
and for different materials of the mechanical element as calculated in the proposed model for gravitational decoherence. In 
addition, Dec(A|_E) p is plotted for the case where gravitational decoherence is not taken into account. When the gap between 
the predicted values with and without gravitational decoherence is large enough, the decoherence estimation formalism allows 
for a test that potentially falsifies the proposed model for gravitational decoherence. The calculations have been made for the 
example experimental parameters go = 1 s _1 , u> m = 1 s _1 and 7 m = 1CP 10 s _1 . 


ence observed - for any no-signalling theory. This means 
that even if quantum mechanics would indeed need to 
be modified we can still draw conclusions from the data 
we obtain. As such, the observations made in such an 
experiment establish a fundamental limit on decoherence 
no matter what the theory might actually look like in 
detail. It is clear, however, that the bound thus ob¬ 
tained is much weaker than if we had assumed quantum 
mechanics. No-signalling is but one of many principles 
obeyed by quantum mechanics, and these other features 
put stronger bounds on the values that Dec(A|.E) can 
take. Our motivation for considering theories which are 
only constrained by no-signalling is to demonstrate even 
such weak demands still allow us to draw meaningful con¬ 
clusions from such an experiment. One can easily adapt 
our approach by introducing further constraints on the 
probabilities Pr[a, 6, c\x, y, z] - but not all of quantum me¬ 
chanics - in order to get stronger bounds. In this case, one 
can similarly obtain an upper bound on Y)ec(A\E) from 
the measured data - this time for the more constrained 
theory. Also in a fully quantum mechanical world, our 
approach yields to a bound (see Figure 3). If we as¬ 
sume quantum mechanics, we may of course also try and 
perform process tomography in order to determine the 
decoherence process, and indeed any experiment should 
try and perform such a tomographic analysis whenever 
possible. The appeal of our approach is rather that we 
can draw conclusions from the experimental data while 
making only very minimal assumptions about the under¬ 
lying physical theory. 


One may wonder, why we only upper bound Dec(A|£’). 
Note that from our experimental statistics we can only 
make statements about the overall decoherence observed 
in the experiment, namely the gravitational decoherence 
(if it exists) as well as any other decoherence introduced 
due to experimental imperfections. Finding that the Bell 
violation is low (and thus maybe Dec(A|.E) might be 
large) can thus not be attributed conclusively to the grav¬ 
itational decoherence process, making a lower bound on 
Dec(A|£() meaningless if our desire is to make statements 
about a particular decoherence process such as gravity. 

Second, we observe that our approach can rule out 
models of gravitational decoherence but not verify a par¬ 
ticular one. It is important to note that a model for 
gravitational decoherence does not stand on its own, but 
is always part of a theory on what states, evolutions and 
measurements behave like. Given such a physical the¬ 
ory and a model for gravitational decoherence, we know 
enough to compute Dec(A|A). In addition, we can com¬ 
pute an upper bound b(j3) on Dec(A|£’) specific to that 
theory, which may give a much stronger bound than no¬ 
signalling alone. Indeed, we see from Figure 3 that this 
is the case for quantum mechanics. Given the calculated 
Dec(A|£’) and the experimentally observed value for b(/3), 
we can then compare: If Dec(A|£’) > b(/3), then the 
model (or indeed theory) we assumed must be wrong. 
However, if Dec(A|.E) ^ 6(/3), then we know that the 
model and theory would be consistent with out experi¬ 
mental observations. We discuss this in more detail in 
the appendix with a candidate decoherence model that 










has been proposed and which - if it is valid - may be 
observed in the experiment suggested above. 

Our approach thus provides a guiding light in the 
search for gravitational decoherence models. It is very 
general, and could in principle be used in conjunction 
with other proposed experimental setups and decoher¬ 
ence models. In particular, it could also be used to 
probe decoherence models conjectured to arise from de¬ 
coherence affecting macroscopic objects, where there ex¬ 
ist proposals to bring such objects into superposition [6]. 
Clearly, however, probing such models using entangle¬ 
ment is extremely challenging. 

It is a very interesting open question to improve our 
analysis and to apply it to other physical theories that 
are more constrained than no-signalling, but yet do not 
quite yield quantum mechanics. Candidates for this may 
come from the study of generalized probabilistic theo¬ 
ries where e.g. [36-41] introduced further constraints in 
order to recover quantum mechanics, but also from sug¬ 


gested ways to modify the Schrodinger equation in order 
to account for non quantum mechanical noise. Since our 
approach could also be applied to higher dimensional sys¬ 
tems, and other Bell inequalities, it is a very interesting 
open question whether other Bell inequalities could be 
used to obtain stronger bounds on Dec(A|.E) from the 
resulting experimental observations. 
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APPENDIX 

Conventions 

For this document, we make the following conventions. 

• The logarithm is with respect to base 2, i.e. log = log 2 . 

• Hilbert spaces are assumed to be finite-dimensional, unless otherwise stated. 

• We denote the set of density operators (states) on a Hilbert space H by H. 

• We identify operators on Hilbert spaces with their reordered versions resulting from permutations of systems. 
For example, for Hilbert spaces Ha, Hb, He and states $ae G S{Ha ® He), &b G S(Hb), we identify the 
state &ae <S> o\b G S(Ha ®He ®Hb) with the state in S{Ha ®Hb®He) resulting from the application of the 
braiding map Ha ®He ® Hb —> Ha ®Hb ®He on ® <tb- 

• For a state pabe € S(Ha ®Hb ®He), we denote its reduced states by according changes of the subscript, e.g. 
Pa '■= Tr b(pae), Pa '■= Tr B ( Tr e {pabe))- 

• For a state pabe G S{Ha ®Hb ®He), entropies are evaluated for the according reduced states, e.g. H(A\B) P 
is the conditional von Neumann entropy of pab = Tr e(Pabe) (c.f. equation (A.8)). 


Appendix A: Background: Decoherence in quantum theory 

In this section, we give a short introduction to decoherence in quantum theory. It consists of concepts, results 
and quantities that are well-established in quantum information science [16]. The topics are chosen to facilitate the 
understanding of our contributions in equations (A.16) and (B.56) rather than to give a full introduction to the 
subject of decoherence. In equation (15), we describe how the dynamical evolution of a system gives rise to a state 
of a tripartite system. This tripartite state plays a central role in our later analysis. In equation (A.8), we explain 
why the min-entropy is the relevant quantity in the information theoretic analysis of decoherence. The min-entropy 
is the quantity that we use for our analysis in equation (A. 16). It is also the quantity that serves as our motivation to 
define a decoherence quantity for generalized probabilistic theories in equation (B.56). We note that a generalization 
of quantum theory by, for example, introducing additional terms into the Schrodinger equation fall under the regime 
of generalized theories in our discussion. 


1. Dynamical evolution and its tripartite purification 

Interaction and non-unitary evolution: Suppose that a system S, initially in a state described by a density 
operator p$ €. S(Hs), undergoes a dynamical evolution over some time interval. If S undergoes this evolution as a 
closed system, then according to one of the postulates of quantum mechanics, the state transforms as 

Ps | -t Us^s Ps Ug^s (A- 1 ) 

for a unitary Ug^-s '■ %s —> Hs (see Figure 6 (a)). In general, however, the system S may be open, i.e. it may interact 
with another system E that is called the environment. We consider the environment E to consist of all the systems 
that interact with system S. Taken together, the combined system SE then forms a closed system and hence evolves 
as 


ps ® Pe Use^se{ps ® Pe)U' se ^ se , (A.2) 

where pE is the initial state of the environment and Use^-se '■ Hs ® He —► Hs ® He is a unitary. 

We may be ignorant about the environment E and only have access to system S. Our description would then treat 
the state of the subsystem S after the evolution as a function of the state ps of S before the evolution. We arrive at 
this description by taking the partial trace over E in expression (A.2): 

Ps | -t Ti\e (Use^se{ps ® PE)Ug E ^ SE j O s^s(ps) ■ 


(A.3) 
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(a) Closed system 



(b) Open system, fixed visible 
part 



(c) Open system, changing 
visible part 


FIG. 6: Different kinds of dynamical evolution, (a) If a system S evolves as a closed system, its state transformation is 
described by a unitary Us-tS- (b) If it evolves in interaction with another system E that is not part of our description, its state 
transformation is described by a TPCPM Qs^s- (c) If the visible part of the overall system before and after the evolution is 
not the same, we describe the state transformation by a TPCPM ©s_>b. 


A map Os-j-s of the form (A.3) is easily shown to be a trace-preserving completely positive map (TPCPM). Thus, the 
evolution of an open system S, when the environment E is not visible, is described by a TPCPM @s-tS ( see Figure 6 

(b))- 

In a yet more general case, it may be that after the evolution of the system SE, we do not have access to system S 
but to a different subsystem B of SE. An example would be a two-particle system S interacting with another two- 
particle system, where we only have access to one particle ( B ) of the four particles after the evolution. Mathematically 
speaking, the fact that we see a different subsystem before and after the evolution means that our factorization of 
the overall Hilbert space changes: Before the evolution, we write H = Hs <S> He' and after the evolution, we write 
H = Hb <8> He ( see Figure 6 (c)). Thus, the unitary evolution of the closed overall system is described by a unitary 
Use'^be '■ Hs <S> He 1 ~> Hb <8> He- Describing only the accessible part before and after the evolution, we end up 
with a TPCPM 

ps Tr B (Use'^,be{ps <S> PE')Ug E ,^ BE ^J =: <3s^b{ps) ■ (A.4) 

Thus, the evolution of a system S to a system B, when being ignorant about the environment, is described by a 
TPCPM Q s ^b- 

Stinespring dilation: We have demonstrated that unitaries on two (or more) systems give rise to TPCPMs on one 
system. It is well-known that the converse is also true: Every TPCPM 0 s_j.b : S(Hs ) —t S(Hb) can be extended 
to a unitary Use'^be on a larger system in the following sense. For Hilbert spaces He' and He of appropriate 
dimensions, it holds that for every pure state |0)(0|s' on He there is a unitary 

Use'^be '■ Hs ® He' —> Hb ® He ( A . 5 ) 

such that 

Tr E {UsE'^BE(ps ® \ty(Q\E')Ug E ,^ BE ) = Q s ^b(ps) for all pg € S(Hs) ■ (A.6) 

This (or an equivalent statement) is the Stinespring dilation theorem [42]. For more details see [16]. 

Textbook definitions of decoherence: From now on, we will take the viewpoint that the TPCPM ©s_j.b is what 
we are given in the first place. Physically speaking, we assume that we are in the setting where all that we observe 
is a process in which a system S in some state ps transforms into a state pb of some system B. We think of this as 
a channel Qs^S: into which we input a system S and get a system B as an output. Our goal in equation (15) is to 
find a precise mathematical formulation of the following question in the quantum theoretical framework: How much 
does the channel decohere the system? 

For the case where S = B, the standard quantum mechanics literature gives some simple descriptions of what the 
decoherence of a system under a dynamical evolution is. As an example, consider the case where A is a spin-1/2 
particle, initially in the spin “up” state in the ^-direction, 

I tz) + I I z) 

V2 


PS = I tx)(tx I , It*) 


(A.7) 
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If the channel Qs^S is given by a measurement of the spin in the 2 -direction, then, written in the 2 -basis, the state 
of the system transforms as 


PS 


1/2 1/2 \ / 1/20 
1/2 1/2 J ^ V 0 1/2 


®s^s(ps) ■ 


(A.8) 


One possible observation one can make in (A.8) is that the spin measurement in the 2 -direction causes the off-diagonal 
terms of the density matrix to vanish. This is an extreme case of the dephasing channel in the 2 -basis, which causes 
a loss of the phase information of the superposition (A.7). This loss of phase information is often equated with 
decoherence. Another feature of (A.8) that is often said to be the characteristic of decoherence is that 0s->s turns 
an initially pure state into a mixed state. 

These descriptions of decoherence, valid in their own right, are not favored by us for mainly three reasons. Firstly, 
these are no quantitative measures of decoherence. Secondly, they lack a clear operational meaning. Thirdly, they 
rely on the quantum mechanical formalism, in which states are expressed as density operators. It is not clear how to 
express them in more general cases that are not described by quantum theory. 

In quantum information science, it is very popular to think of the systems arising in the purified picture we just 
presented as being controlled by parties with intentions and interests rather than just being dead physical objects. 
We will follow this spirit and from now on use the language of a game and speak of parties Alice, Bob and Eve, that 
we think of as agents controlling the systems A, B and E. 


2. The min-entropy as a measure for decoherence 

The coherent information: As mentioned in equation (15), it has been realized in quantum information science 
that important quantitative measures of the channel are functions of the state pabe that we described above. One 
such measure quantifying decoherence is the coherent information [43]. It is defined in terms of the conditional von 
Neumann entropy 


B(A\B) P :=B(AB) P -B(B) P , (A.9) 

where H(AB) p = — Tr(pAB log {pab)) and H(B) = —Tr(p B log(p B )) is the von Neumann entropy of the reduced state 
Pab and p B , respectively. The coherent information is defined as 

I(A)B) P :=-B{A\B) P . (A.10) 

The coherent information I{A)B) p has been shown to be related to the quantum channel capacity Q(Qs^b) of O 
which is known as the Lloyd-Shor-Devetak (LSD) theorem [44-46]. It says that 

Q(©s^b) = lim - max I(A n )B n ) p , (A.ll) 

n-s-oo n p s neS(H.f n ) 

where I{A n )B n ) p is the coherent information for pA n B n = l® ra <S> 0f// B (PA n s n ) and pa^S 71 is a purification of psn. 
The state ® Of/ >B (p J 4 " 5 ") results from the n-fold use of the channel O B _> B to transmit S n , i.e. n copies of 
system S, while the purification A n of S n remains unchanged. Thus, the r.h.s. of (A.ll) is the coherent information 
in the limit of infinitely many channel uses. Likewise, the quantum capacity <3(0 s->b) is the limit of the achievable 
rate for quantum data transmission in the limit of infinitely many channel uses. One says that the quantum capacity, 
and therefore the coherent information, is an asymptotic quantity. This has the disadvantage that from the coherent 
information, only very limited statements about finitely many uses of the channel can be made. 

The min-entropy: More insight about the behavior of the channel under finitely many uses can be gained by 
considering the corresponding single-shot quantity. To formulate it, note that the state pabe is pure, in which case 
the duality relation H{A\B) p = — B{A\E) p for the conditional von Neumann entropy holds. This gives us 

I{A)B) P = B{A\E) P . (A. 12) 

The corresponding single-shot quantity for the conditional von Neumann entropy B{A\E) p is the conditional min- 
entropy, or just min-entropy, B m [ n (A\E) p [20]. It is defined as 

B min (A\E) p = max sup {A G R. | pae 2~ x t A <g> cr E } , 

<7 E 


(A. 13) 
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where the maximum is taken over all subnormalized density operators on He, he. all positive operators on He with 
trace between 0 and 1. The min-entropy quantifies the maximal size of a subsystem of A that can be decoupled from 
E [17], and thus tells us how many EPR pairs between Alice and Bob can be created [47] given a noisy output state 
Pab- To obtain the single-shot capacity of n channel uses we are - as in the asymptotic case - allowed to optimze over 
input states pA n S n - Clearly, however, the resulting expression can be lower bounded using a particlar input state given 
by n copies of the maximally entangled state. This is the test state we employ here, and hence our test also provides 
a bound on the single shot capacity. For instance if A is a 2 level system, then the min-entropy readily quantifies the 
number of EPR-pairs we can recover, given that we started with n EPR pairs as an input. The min-entropy thus has 
a very appealing operational interpretation. 

For our purposes, another expression for the min-entropy is more useful. In the following, we use the symbol ~ to 
denote that two Hilbert spaces are isomorphic, i.e. Ha — Ha> means that the two spaces have the same dimension. 
It has been shown [19] that the min-entropy can be expressed as 

H min (A\E) p = —logoU max F 2 ($ A A', 1a ® TIe^A'(pae)) , (A. 14) 


where d>A is the dimension of the Hilbert space Ha of system A, A' is a system with H a' — Ha, the maximization is 
carried out over all TPCPMs 7 Ze^a' from system E to system A', F(p, a) = Tr \Jp l / 2 ap v / 2 is the fidelity and $, 4 ^ 4 / 
is a maximally entangled state on AA' , i.e. <1 >aa' is an element of the set 


Taa' 


| $4A' € S (Ha <S> 'Ha' ) 


There are bases {|i)^}j, {| i)A'}i of Ha, Ha 1 such that <&aa' = 
\4>){4>\aa' with | 0 >aa' = ^2 ® 


(A. 15) 


The choice of <&aa' €E Taa', i.e. the choice of bases for Ha and Ha', is irrelevant for the value of H m i n (A\E) p . Since 
every $, 4 , 4 / = |<^)(^>|aa' € Taa' is pure, we have that F(<S>aa', &AA') = \/(</>| <t\4>)aa' for any state ctaa' on AA'. 

The expression (A.14) provides an intuition for the min-entropy. We think of the system ABE, which is in the 
pure state pabe , as being distributed between Alice, Bob and Eve. Imagine that Eve tries to perform operations on 
her share of the system with the intention to bring the reduced state between her and Alice as close as possible to 
the maximally entangled state $, 4 , 4 /, where the square of the fidelity is the measure of closeness. The closer Eve can 
bring the state to the maximally entangled state, the smaller the min-entropy H m i n (A\E) p . The overall situation of 
our decoherence analysis is shown in Figure 7. 


E' 




A' 


B 

A 


FIG. 7: Overall picture of our decoherence analysis for quantum theory. We interpret H mln (A\E) p as a measure for 
how close Eve can get to maximal entanglement with Alice by applying a transformation 7 Ze->a' 011 her share of the state 
Pabe- 


The min-entropy is strictly more informative than the conditional von Neumann entropy in the following sense. In 
the iid limit (which stands for independent and identically distributed), where many identically prepared systems go 
through the channel and end up in a state the min-entropy converges to the conditional von Neumann entropy: 

lim -HF n (A n \B n ) p ® n = H(A\B) p , (A.16) 

n—>-oo 71 

where e > 0 is an arbtirary smoothing parameter. This is known as the asymptotic equipartition property [48]. Thus, 
in the limit of infinitely many channel uses, where the asymptotic quantity is relevant, the min-entropy reproduces 
the conditional von Neumann entropy. 

To gain some intuition for H m i n (A\E) p , we now have a look at some special cases. For these special cases, we 
assume that Hs — Ha — Hb — He- Assume that initially, the state pas is maximally entangled, i.e. pas = 4>a.s for 
some $^45 £ F^g analogous to (A.15). We think of the channel purification Use'^be as being controlled by Eve. 

• If the adversary Eve leaves system S untouched, i.e. the channel Qs^b is the identity channel (or any other 
unitary channel), then pae = Cl^a <8> Pe for some state pE of system E. In that case, H min (A\E) p = log <7,4, 
and we say that there is no decoherence. 
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• In the other extreme case, Eve snatches away the system S and forwards an uncorrelated system to Bob. In 
this case, pae £ r ae with analogous to equation (A. 15) (maximal entanglement between A and E ). Then, 
Hmin(A\E) p = — logd^, and we say that we have full decoherence. 

• As an intermediate case, we might consider the case where Eve interferes such that she does not end up with 
maximal entanglement with Alice but such that she is classically correlated with Alice in some basis, i.e. 
Pae = d-A 1 J2k \ k )( k U ® \k)(k\ E - In that case, H min (A\E) p = 0, and we speak of partial decoherence. 


Appendix B: Decoherence estimation through CHSH tests in quantum theory 


1. Introduction 


Our goal is to show that Alice and Bob can estimate the decoherence by performing a Bell experiment. We pose 
it as a feasibility problem: is it possible to observe certain statistics in a Bell experiment given a certain level of 
decoherencel Solving this problem allows us to determine and plot the feasible region in the space of suitably chosen 
parameters. 

We look at the simplest Bell experiment, known as the Clauser-Horne-Shimony-Holt (CHSH) [33] scenario. If pab 
is the state that Alice and Bob share and Aj,B% for j,k £ {0,1} are the observables they perform, then the CHSH 
value equals 

/3 = Tr([A 0 <8) B 0 + A 0 ® B 1 + A x <g> B 0 — Ar <g> B^pab). (B.l) 

As explained previously the min-entropy H m i n (A\E) defined in Eq. (A.13) captures the notion of decoherence 
between Alice and Bob (although note that high min-entropy corresponds to low decoherence and vice versa). Since 
the range of values that the min-entropy takes depends on the dimension of Alice’s system (denoted by d^), it is only 
meaningful to compare scenarios in which is fixed. For simplicity, we consider the simplest non-trivial scenario in 
which the subsystems held by Alice and Bob are qubits, dA = ds = 2. 

We define the feasible region S as follows. A pair of real numbers (u, v), where u £ [—1,1] and v £ [0, 2-\/2] belongs 
to S if there exists a tripartite state pabe and binary observables Ao,A\ on Ha and Bq,Bi on Hb such that 

• subsystems A and B are qubits: dirndl = dim'Hs = 2 

• The conditional min-entropy of A given E equals u: H m - ln (A\E) = u. 

• The CHSH value given by Eq. (B.l) equals v: f3 = v. 

First note that a CHSH value of v if 2 can be achieved using trivial measurements (namely {1,0}) acting on an 
arbitrary state. Therefore, for v ^ 2 all values of u £ [—1,1] are allowed. For the remainder of the argument we 
implicitly assume that v > 2 and the following intuitive argument shows why certain pairs (it, v) must indeed be 
forbidden. Consider a point u ss — 1 and v > 2. According to the operational meaning of the min-entropy (A. 14), 
it ss — 1 means that Eve can recover the maximally entangled state with Alice with fidelity close to unity, which 
clearly allows Alice and Eve to violate the CHSH inequality. On the other hand, since v > 2 Alice also observes a 
CHSH violation with Bob. This violates the monogamy relation for tripartite three-qubit states proved in Ref. [49], 
which states that Alice can violate the CHSH inequality with at most one party (even if she is allowed to use different 
measurements for different scenarios). This simple argument leads to the conclusion that the region u « — 1 and v > 2 
is forbidden. In the remainder of this section we show that the non-trivial part of the feasible region S can be fully 
characterised by a single inequality. 


Theorem B.l: A pair of real numbers (it, u) where u £ [—1,1] and v £ (2,2y/2] belongs to the feasible region S if 
and only if 


where 


u > /( v), 

f{v) := 3 - 2 log max ^ 2a/1 + c z + ^ 1 - c z + -^= 


where the maximization is taken over 


v 

7T 



(B.2) 


(B.3) 


-1 < c, < 1 - 


(B.4) 
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While the definition of / might seem complicated, it is straightforward to see that / is monotonically increasing in v 
and evaluating f(v) numerically for a particular value of v is straightforward since the function to be maximized is 
concave. The feasible region S is plotted below. 

The proof of Theorem B.l is conceptually simple, but it requires a wide array of technical tools, which we present in 
equation (B.4). In equations (B.52) and (B.54) we prove the direct and converse parts of Theorem B.l, respectively. 


2. Preliminaries 

Definition B.2: Let Ha, Hb be Hilbert spaces of dimension d. A generalized Bell basis for Ha ®Hb is a set 
{1}^_i of d 2 pure states on Ha ®Hb which satisfy 

Tr A \$ j )($ j \ = ^, Tr B |^)(^| = ^ for j = 1,..., d 2 and (B.5) 

d 2 

| = 1 A 0 1 B . (B.6) 

A state pab G S(Ha ®Hb) is called Bell-diagonal if it is diagonal in some generalized Bell basis, i.e. if there exists 
a probability distribution {pj}y =1 such that 

d 2 

pAB = Y,Pf l*iX*il- ( B - 7 ) 

1=1 

Lemma B.3: Let A ^ 0 be a positive semi-definite operator , let II a be the projector on its support and let |6) be a 
normalised vector. Then A ^ | b) (b \ iff 

n A \b) = \b) and (blA^lb) < 1. (B.8) 

Note that since A might not be invertible, A _1 is only defined on the support of A. 


a. Two-qubit states 


A two-qubit state written in the Pauli basis takes the form 

1 
4 


Pab = j{^a 0 Is + ^2 a j a j 0 Is + ^a 0 + '^ J T jk°j 0 &k), 

i 1 l,fe 


(B.9) 




16 


where all the summations go over {x, y, z}. It is known that for every state there exists a local unitary Ua <8> U B 
which diagonalizes the correlation tensor (i.e. ensures that Tjk = 0 for j ^ k) and since all the properties we consider 
are invariant under local unitaries we can make this assumption without loss of generality. We denote these diagonal 
entries T xx , T yy and T zz by c x ,c y and c z , respectively, which simplifies the expression to 

pAB = “ (H-A (8) 1b + ^2/CLjGj 0 t B + &A ® ^ bjCTj + ^2cj(Tj 0 Gj). (B.10) 

3 3 3 

Without loss of generality, we assume that | c x \ ^ \c y \ ^ \c z \ and c x ,c y ^ 0. As shown in Ref. [50] every Bell-diagonal 
state of two qubits (up to local unitaries which, again, we can safely ignore) can be written as 


PAB = 


(B.ll) 


1=1 


where {pj}q =1 is a probability distribution and |t&i 2 ) = anc l |^ 3 , 4 ) = ^ 01 ^ 10 ^ ■ It is easy to verify that 


\/2 


Pab = - (1a <8> Is + 


CjVj <8 <Tj), 


(B.12) 


where 


c x = Pi -P 2 +P 3 -P 4 , 

Cy = ~Pl +P2+P3-P4, (B. 13) 

C z = Pi + P2 ~ P3 - P4 ■ 


b. Non-locality 

Definition B.4: For a bipartite quantum state pab the maximum CHSH value is defined as 

Pm ax (PAB) ■= A max Tr[(A 0 ®B 0 + A< ) <8)B 1 +A 1 ®B 0 -Ai®B 1 )p AB ], (B.14) 

A 0 ,A 

where the maximisation is taken over all Hermitian, binary observables. 

Note that for all states /3 max ^ 2 and we say that the state violates the CHSH inequality if /3 max > 2. It was shown 
in Ref. [51] that if pab is a state of two qubits then the value of /? m ax is fully determined by the correlation tensor. 
Adopting the convention |c x | ^ |c y | ^ |c 0 | we have 


Ana x(PAb) = 



if c£ + c£< 1, 
otherwise. 


(B.15) 


c. Entropic measures of entanglement 

To derive a bound on the min-entropy H m i n (A\E) p , we will use a closely related quantity, namely the max-entropy. 

Definition B.5: For a bipartite quantum state pab the conditional max-entropy (or just max-entropy) is 
defined as 

H max (A\B) = max log d A F 2 (pab , ^A ® cr B ), (B.16) 

0’S 

where it a is the maximally mixed state on A and the maximisation is taken over all states on B. 

The proof uses the following known properties of the min- and max-entropies. 
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Lemma B.6 (Duality. [19]): Let pabc be a tripartite state. Then 

H max (A\B) p + H min (A\C) p >0, 

and the equality holds iff pabc is pure. 

Lemma B.7 (Data-processing inequality, [20]): For an arbitrary tripartite state pabc we have 

H ma JA\B) 3s H max {A\BC). (B.17) 

Lemma B.8 (Conditioning on classical information, Proposition 4.6 of [52]): Let pabi< be a tripartite state where 
K is a classical register: 

PABK = E VkT k AB ®\k){k\. (B.18) 

k 

Then 

H max (A\BK) PABK = log (^p fc 2 ffm “ (A|B) ^). (B.19) 

k 

Finally, we need an explicit expression for the max-entropy of a Bell-diagonal state. Note that by assumption 
dA = ds = d. 

Lemma B.9: Let pab be a Bell-diagonal state of form (B.7). Then the conditional max-entropy equals 

J ffmax(4l|B) = -log(i + 21og(^ v ^-). (B.20) 

3 

To prove Lemma B.9 we use the fact that the optimization problem which appears in the definition of the max-entropy 
(B.16) can be written as a semidefhhte program (SDP) [53]. More specifically, given pab we have Lf max (A\B) = log A, 
where A is the value of the following SDP for pabc being an arbitrary purification of pab 


PRIMAL : minimize p 

DUAL : maximize L^{pabcYabc) 

subject to pts ^ tr a(Zab) 

subject to Tr c(Yabc) ^ 1 a®c B 

Zab <8> 1 c ^ pabc 

Trcr B ^ 1 

Zab € V(LLab) 

Yabc € VILLabc) 

p ^ 0 

<jb € V(LLb) 


where 'P('H) denotes the set of positive semi-definite operators acting on LL. By providing feasible solutions for the 
PRIMAL and the DUAL we show that for Bell-diagonal states 

A = ^(EvE) 2 ( R21 ) 

3 

which is precisely the statement of Lemma B.9. 

Proof. Let pabc = \'f>ABc){4>ABc\ be a purification of pab, e.g. 

\iPabc) = E v'/'/ <l b) ® b’>- (B.22) 

3 

For the PRIMAL consider 

Zab = (£ -y/Pfcl ^fc)(^fc |) (B.23) 

3 


(B.24) 
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Clearly, Zab ^ 0, p ^ 0 and since Tr a{Zab) = (Xq the constraint is easy to check. The last 

inequality we need to check is 


(Ev^) E v / Pfcl <J) fc)( <J> fcl ® Ac > Pa.bc- (B.25) 

j fc 


We apply Lemma B.3 to A = Zab <S> A c and \b) 


I i>ABc)- The projector on the support of Zab ® A c equals 


n= E I *,■><*,• I ®ic 

j'-Pj >o 


(B.26) 


and it is easy to verify that II| rpABc) = IVabc)- Moreover, since ( Zabc ) 1 = {Zab) 1 ( g > Ac we have 

=(E«r e VPm(®m\ ® (m\ E i $ ») ® h 

i \ m ) \n:p n >0 / 


= 1. 


Showing that .Zab and p constitute a valid solution to the PRIMAL implies that A ^ g( ]TT y/PjY 
For the DUAL consider 


Y ABC = ±Y,\*M**\®\m, 


jk 


<J B = 


11 B 

~d' 


Note that Yabc is proportional to a rank-1 projector. The first constraint gives 

Trc(TABc) = ~j E I I = ^ 1a <8 1b = 1a <8 cr B 

j 


(B.27) 

(B.28) 


(B.29) 

(B.30) 

(B.31) 


and the remaining ones are easily verified to be true. The value of this solution equals Tr {pabcYabc) = g ( Xq y/Pj) 
which implies that A ^ ^ y/Pj ) 2 ■ □ 


d. Sufficiency of considering Bell-diagonal states 

To prove the converse part of Theorem B.l, we will use the following argument, which is similar in spirit and inspired 
by the symmetrization argument presented in Ref. [23]. 

Lemma B.10: Let pab be an arbitrary state of two qubits. Then, there exists a Bell-diagonal state oab which 
satisfies 

/3ma x(pab) = /?ma x(cab) and H max (A\B) a > H max (A\B) p . (B.32) 

Proof. We present an explicit construction of <tab which meets the requirements. According to Eq. (B.10), pab can 
be written as 


Pab = ^ (l a <8 1s + E “TT? ® As + A a 8 E ? b cr j + E c i <J i ® a j) • 
j j j 

Moreover, consider the following random unitary channel 

1 4 

A(pab) = - YPi ® U Ypab{U} ® U}), 
i =1 


(B.33) 


(B.34) 
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where U\ = 1, U 2 = g x , U 3 = a y and U 4 = a z . It is easy to verify that for j £ {x, y , z} 


A ((Tj 0 Is) = A(l yi 0 aj) = 0 (B.35) 

because each Pauli operator commutes with identity and itself but anticommutes with the other two unitaries. This 
implies that gab = A-(pab) is Bell-diagonal. Moreover, one can check that the map preserves the correlation tensor, 
i.e. for j £ {x, y, z} 


A(Gj 0 Gj) = Gj 0 Gj, (B.36) 

which implies that P max (pAB) = Pmax(o'AB)- To check the last property consider the following state 

1 4 

&abk = ^ S> Uj)pA B (Uj 0 U}) 0 \j)(j\. (B.37) 

1=1 

By the data processing inequality, we have H max (A\B) rJ ^ H max (A\BK) a and by conditioning on classical information 
we have 

4 1 

H max (A\BI<) a = log ( J] 4 • ) , (B.38) 

4=1 4 

where r j AB = (Uj 0 Uj)pAB(Uj 0 Uj). (B.39) 

Since the max-entropy is invariant under local unitaries we have H max (A\B) r j = H max (A\B) p for j £ {x,y, z} which 
implies that 


H max (A\B) a ^ H max (A\BK) a = H max (A\B) p . (B.40) 

□ 

The final technical lemma concerns the problem of maximizing the max-entropy of a Bell-diagonal state of two qubits 
whose maximal CHSH violation is fixed. 

Lemma B.ll: Let pab be a Bell-diagonal state of two qubits, whose maximal CHSH violation equals (3 £ (2,2y/2\. 
Then, the max-entropy of pab satisfies the following inequality 

H max {A\B) < -f(P) (B.41) 

for function f defined in Eq. (B.3). Moreoever, there exists a state which saturates this inequality. 

Proof. According to Lemma B.9 the max-entropy of a Bell-diagonal state of two qubits equals 

4 

H max (A\B) = -l + 2log(22VPj)- (B.42) 

l=i 

Here, it is convenient to express the probabilities through the correlation coefficients c x ,c y ,c z . Inverting Eqs. (B.13) 
gives 


Pi = ^(1 + C x - Cy + C z ), 

P2 = ^(1 - C x + Cy + C z ), 

(B.43) 

P 3 = ^(1 +C X 1 Cy - c z ), 

Pa = ^ ( 1 c x Cy c z ), 

(B.44) 


which allows us to write 

H max (A\B) = -3+ 2 log g{c x ,c v ,c z ), (B.45) 


g(c x ,c y ,c z ) — \J\ + c x — c y + c z + \J\ — c x + c y + c z + \J\ + c x + c y — C z + \J\ — c x — c y — C z . (B.46) 


where 
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In the space of correlation coefficients the feasible set are the triples ( c x ,c y ,c z ) for which the function g(c x ,c y ,c z ) is 
well-defined (the expressions under the roots must be non-negative). As before, we assume without loss of generality 
that \c x \ ^ \c y \ ^ \c z \ and c x ,c y ^ 0. Then, the maximal CHSH violation (we are only interested in states that violate 
the CHSH inequality) is given by Eq. (B.15) 


P = 2 \J c l + c l' 

Since in our case j3 is fixed, the angular parametrisation takes the form 

Cx = ~sp2, ^ + D and ° v = 71 C ° S + D’ 

where q = ^ and cj> £ [0, 7t/ 4] (which ensures c x ^ c y ^ 0). Note that 

c x +Cy = q cos cj>, 

Cx- c y = q sin cj>. 

It is easy to check that the allowed range of c z is 

q sin cj) — 1 ^ c z ^ 1 — q cos cj>. 

Note that we should also impose the condition \c z \ ^ \c y \ but as it turns out the optimal solution will satisfy it even 
if we do not include it explicitly. To maximize the max-entropy it is sufficient to maximize function g defined in 
Eq. (B.46), which in the angular parametrisation equals 

g(cj>, c z ) = y/l + c z + qsincj) + y/l + c z - qsincj) + yA - c z T qcoscj) + yjl - c z - qcos(f>, (B.47) 


over 


7Z = {(</), c z ) : <j> G [0, 7 t/4 ], q sin cj) — 1 < c z < 1 — qcoscj)}. 


(B.48) 


The maximum is achieved either in the interior (denoted by TZ ln t) or at the boundary. Let us start by ruling out the 
first option. Function g is differentiable everywhere in lZi n t and the partial derivatives are 


dq 1 1 -1 -1 

dc z 2 -yH + c z + q sin cj) 2-^1 + c z — qsincj) 2 y/T^c^+~q~cos~cj) 2^/1 — c z — qcoscj 

dg qcoscj) —qcoscj) —qsincj qsincj 

dcj) 2 -\/l + c z + qsincj) 2 yJ\ + c z — qsincj 2 yJ\ — c z + qcoscj) 2-^1 — c z — qcoscj) 


(B.49) 

(B.50) 


To prove that there is no maximum in the interior, it suffices to show that there is no ( cf>,c z ) G TZ m t such that both 
derivatives vanish = 0. To do this we consider the following linear combination 

.. . „ . . dg 2 dg sin cj) + cos <f> sin cj — cos cj —2 sin cj 

s(4>, c z ) = 2 sm cj> ■ — -1— • — = -|— . H- . 

dc z q ocj y/1 + c z + qsincj) yl + c z — qsincj) yl — c z + qcoscj) 


and show that s(cj,c z ) = 0 has no solution in 7Z ln t- Since the last term of s(cj,c z ) is negative, a necessary condition 
for s(cj, c z ) = 0 is that the sum of the first two terms is non-negative, which is equivalent to 


sm < 


cost 




cos o — sm t 


y /1 + c z + qsincj> ^ ^1 + c z — <7 sin 


This can be rearranged to give 


c z > 


q 


- 1 , 


2 cos cj) 

which contradicts the second inequality in the definition of 7Z m t as shown below. 

c z ^ —-—- — 1 and 1 — q cos cf> > c z 
2 cos cp 

q 12 

=4> 1 — q cos cj) > - -- - 1 --- + cos </><-. 

2 cos cj) 2 cos cj) q 


(B.51) 

(B.52) 
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It is easy to check that the left-hand side of the final inequality is always at least \/2, while the right hand side is 
always at most y/2. This proves that the final (strict) inequality is always false, which implies that s(</>, c z ) = 0 has 
no solutions in lZi nt and that g(<j>,c z ) has no maximum in lZ m t- 

The boundaries c z = q sin </> — 1 and c z = 1 — q cos <j> correspond to one of the expression under the roots being zero. 
Since the square root function has infinite slope at 0, such solutions cannot be optimal. Therefore, the maximum 
must be achieved at the boundary <fr = 0. Combining Equations (B.45) and (B.47) and setting 0 = 0 leads directly to 
the statement of the lemma. 

To show that the solution of the optimization problem satisfies \c z \ ^ \c y |, it is sufficient to show that for 0 = 0 
and c z = — c y = —q/2 the partial derivative dg/dc z is strictly positive. □ 


3. The direct part 

Here, we show (by an explicit construction) that points described by v € (2, 2-\/2] and f(v) ^ u ^ 1 are allowed. 
Lemma B.ll shows that for v € (2,2\/2] there exists a Bell-diagonal state of two qubits whose max-entropy equals 

H m UA\B) = -/(«). (B.53) 

By duality (Lemma B.6), if pabe is an arbitrary purification, the conditional min-entropy equals 

H min (A\E) = f(v). (B.54) 

In this example u = f(v), which corresponds to a point lying precisely on the boundary defined in Theorem B.l. In 
order to obtain higher values of u (all the way up to 1), it suffices to apply noise of appropriate strength to subsystem 
E. 


4. The converse part 

Here, we show that every feasible point (u. v) must satisfy u ^ f(v). Consider a state pabe for which H m i n (A\E) p = u 
and which for some measurements achieves the CHSH value of v. Clearly, /3 ma x(pab) v and by Lemma B.6 
H max (A\B) p ^ —u. Applying the symmetrization argument (Lemma B.10) gives rise to a Bell-diagonal state (tab 
such that H max (A\B) a ^ —u and /3 ma x(cab) ^ v. By Lemma B.ll these quantities must satisfy 

Hmax(A|£) CT < -f(p max (<T AB )), (B.55) 

which implies that 


U Si -H max (A\B) a ^ f( 0max{<T Ab)) > f(v), 
where the last inequality follows from the fact that / is monotonically increasing. 


(B.56) 


Appendix C: Decoherence estimation through CHSH tests in GPTs 

In this section, we are going to develop a framework for decoherence analysis in analogy to equation (15), but 
without assuming that nature is correctly described by quantum theory. Instead, we will work in a framework that 
makes only minimal assumptions about the probabilistic structure of measurements. This allows to make statements 
in cases where quantum theory might not be a correct description of nature. 

In equation (B.56), we define a framework for probabilistic theories that has become a standard one in the literature. 
Besides defining the core structure in equation (B.56), we explain in equation (C.8) how we extend this framework 
to make it suitable for analyzing tripartite states, in a way that allows us to make a decoherence analysis that is 
analogous to the quantum case. 

In equation (C.19), we will define a decoherence quantity Dec(A|E') UJ for GPTs as an analogue of the quantum 
min-entropy H m i n (A\E) p . This will be our quantity of interest for the decoherence analysis for GPTs. We will first 
motivate an expression for Dec(A|ii) w in equation (C.19), inspired by expression (A.14) for the min-entropy in the 
quantum case. This expression will require us to say what a maximally entangled state in a GPT is. We will define 
it in equation (C.22). 

equation (C.34) is devoted to finding a bound on our decoherence quantity in terms of the CHSH winning probability 
for Alice and Bob. This is a measurable quantity in the case where the channel is an iid (for independent and 
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identically distributed) channel, meaning that it behaves identically in repeated uses of the channel without building 
up correlations amongst systems going through the channel in different uses of it. This is a practically relevant case, 
giving our bound a practical meaning. This bound allows us to infer non-trivial statements about decoherence from 
measured data when, apart from the iid assumption, we assume only very little about the behavior of nature. We 
approach our bound by first bounding our fidelity-based decoherence quantity by a trace distance-based quantity. We 
will then bound this trace distance-based quantity in terms of the CHSH winning probability for Alice and Bob by a 
quantity that can be expressed as a linear program. 

Finally, in equation (C.89), we show how our bound can be expressed as a linear program and present the numerical 
results. This is followed by a discussion of the physical interpretation of our numerical findings. 


1. The framework 


a. A basic framework for GPTs 


Frameworks for probabilistic theories in which quantum theory and classical theory can be formulated as special 
cases have already been considered some decades ago [54-56]. After some period of oblivion, a seminal paper by Hardy 
[57] caused a revival in the interest in such frameworks (see, for example, [36-41] and references therein). Today, they 
are generally refered to as frameworks for generalized probabilistic theories [58]. 

We formalize our decoherence analysis for GPTs in the abstract state space framework [59-62]. It is one rigorous 
formalization of what a generalized probabilistic theory is, amongst a few equivalent or closely related ones that can 
be found in the literature (see the references cited above). We prefer it for its concise and precise formulation. For 
the sake of brevity, we will not go far beyond the mere mathematical definitions related to abstract state spaces here. 
For a detailed introduction to abstract state spaces, see [63]. 

Definition C.l: An abstract state space is a triple (V, V + 1 u), where V is a finite- dimensional real vector space, 
V + is a cone 2 in V which is closed 3 and generating 4 and u € V* is a linear functional 5 on V such that u(oj) > 0 for 
all oj £ V + \ {0}. The functional u is called the unit effect. 

Definition C.2: For an abstract state space (V, V + ,u), we define the following induced structure (see Figure 8): 
The normalized states are the elements of the set 


!! := {w £ V + | u(oj ) = 1} . 

The subnormalized states are the elements of the set 

:= {u> £ V + | u(lo) ^ 1} . 

The effects are the elements of the set 

£:= {e £ V* | 0 < e(w) < 1 Vw € ft} . 
The measurements are the elements of the set 

A4 := f. M C £ finite 



(C.l) 


(C.2) 

(C.3) 


(C.4) 


An effect respresents a measurement outcome. If a system in a state w is measured with respect to a measurement 
M = {ei,..., e„}, then efc(w) is the probability that the measurement yields the outcome associated with e&. 


2 A subset V - *" C V is a cone in V if 

(Cl) y+ + c v+, 

(C2) aV+ C V+ for all c* ^ 0, 

(C3) fl ( -V + ) = {0}, 

3 We assume the standard topology on V, i.e. the only linear Hausdorff topology on V. 

4 A cone V+ C V is generating if V + — V + = V. 

5 For a finite-dimensional vector space V, we denote by V* the dual space of V , i.e. the vector space of linear functionals on V. 
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FIG. 8: Visualization of an abstract state space. The set V + is a cone in a vector space V, here shown as a cone with a 
square base. The normalized states Q. are given by the intersection of V + with the plane on which the functional u takes the 
value 1. The set fP consists of those elements of V + on which the functional u takes values between 0 and 1. 


Example C.3 (Quantum theory): The probabilistic structure of measurements on a (finite-dimensional) quantum 
system can be formulated as an abstract state space. For a quantum system with an associated Hilbert space H, 
consider the abstract state space (V,V + ,u) = (Herm('H), Pos("H), Tr), where V = Herm(7Q is the real vector space 
of Hermitian operators on H, V + = Pos("H) is the cone of positive operators on H and u = Tr is the trace on H. 
According to Definition C.2, this yields the states fl = {p £ Pos (H) | Tr (p) = 1}, which are precisely the density 
operators on H. Analogously, are the subnormalized density operators. The effects are the functionals induced 
by POVM elements via the trace, £ = {Tr(P • ) | P £ Pos (H),P ^ In}- Accordingly, the measurements are the 
sets of functionals that are induced by POVMs, Ai = {{Tr(P^ ■ ) | P k £ {Pfcjfc} | {Pk}k Is a POVM}. This precisely 
reproduces the structure of measurement statistics in quantum theory. For further details, see [63], ■ 

By our definition, £ is the set of all linear functionals e such that 0 ^ e(w) ^ 1 for all oj £ ft. The underlying 
assumption that every such linear functional represents a physical measurement outcome has been called the no¬ 
restriction hypothesis [40]. A priori, there seems to be no immediate physical reason for this assumption, and some 
authors have argued about how to weaken this assumption [64]. For our purposes here, it is not relevant whether the 
no-restriction hypothesis holds, and weakening the assumption complicates the definitions. Thus, we assume it for 
simplicity. 

In equation (C.19), we will define a decoherence quantity Dec^lP^ analogous to the quantum min-entropy 
^min(A|P) p . We will take our inspiration from expression (A.14) for the quantum min-entropy, which involves the 
fidelity as a measure of closeness of quantum states. Therefore, it is desirable to have a generalization of the fidelity 
to states in abstract state spaces. Such a generalization is easily found once it is noticed that the quantum fidelity of 
two states can be expressed as the Bhattacharyya coefficient (or classical fidelity) of the two probability distributions 
that the two states induce, minimized over all measurements. More precisely, the quantum fidelity satisfies [25] 

F(p, cr) = min E v / Tr(P fc p)Tr(P fc cr), (C.5) 

{ Pk }k 

where the minimization runs over all POVMs {Pk}k on the Hilbert space on which p and a are defined. The sum in 
(C.5) is precisely the Bhattacharyya coefficient of the probability distributions that the POVM {Pfcjfc induces on the 
states p and cr. This motivates us to define the fidelity for abstract state spaces as follows. 

Definition C.4: Let (V, V + ,u) be an abstract state space with normalized states fl and measurements A4. For 
states uj,t £ fl, we define the fidelity of uj and r as 

F(cv,t) := inf b(u),r\M), where 6 (o;,t|M)=V^ y / e(w)y / e(r). (C.6) 

MeM z ' 

eSM 

The quantity b(u, t\M) is the Bhattacharyya coefficient (or sometimes called the classical fidelity) of the probability 
distributions that the measurement M induces on the states u> and r. 

The fidelity as defined in Definition C.4 precisely reduces to the quantum fidelity in the case where the abstract state 
space is a quantum state space. In addition to the fidelity, in equation (C.34) we will also consider a generalization of 
the quantum trace distance D(p, a) = fTr\p— cr| in order to formulate a bound on Dec(A|£') w . Somewhat analogously 
to the fidelity, the quantum trace distance is equal to the total variation distance (or classical trace distance) between 
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the two probability distributions that the two states induce, maximized over all measurements [25]: 

D(p,<r) = max \ |Tr(P fc p) - Tr(P k a)\. 

{.Hkfk ^ , 


(C.7) 


This motivates the following definition. 

Definition C.5: Let ( V,V + ,u ) be an abstract state space with normalized states fi and measurements A4. For 
states w, r £ ft, the trace distance between oj and r is given by 

D(ui,t) := sup d(w,r|M), where d(oj, t\M) — - |e(w) — e(r)| (C.8) 

2 ttu 

The quantity d(w,r|M) is the total variation distance (or sometimes called the classical trace distance) between 
the probability distributions that the measurement M induces on the states lo and r. 

Note that the fidelity and the trace distance take values between 0 and 1 for all states. For squares of the quantities 
F, b, D and d, we will write the square sign right after the letter, e.g. we will write F 2 (w,r) instead of (F(uj,t)) 2 . 


b. A tripartite framework for GPTs 

In equation (C.19), we will consider a tripartite situation for the decoherence analysis, analogous to equation (15). 
This requires us to model a tripartite scenario mathematically since such a structure is not induced by an abstract 
state space (V. V + ,u) alone. We need to specify it as additional structure. Our goal here is to do this with the weakest 
possible assumptions, resulting in a very general validity of the bounds we derive. 

Instead of assuming individual state spaces for every party, we only consider their overall combined state space, 
modelled by an abstract state space (V, V + ,u) and all its induced structure as in Definitions C.l and C.2. This has 
the advantage that we do not have to make assumptions about how individual state spaces combine to multipartite 
state spaces, keeping our assumptions weak. For our purposes, the only structure that we need to add to an abstract 
state space (V,V + ,u) to make it suitable for the description of a tripartite scenario are the local transformations 
that each individual party can perform. The local measurements of the three parties are then induced by these local 
transformations. 

We consider three parties, which we call Alice (A), Bob ( B ) and Eve ( E ) as before. We begin our considerations by 
assuming that there are three sets Ta, Tb and Te, containing all the transformations that Alice, Bob and Eve can 
perform, respectively. By a transformation, we mean a linear map T : V — > V which maps states to subnormalized 
states, i.e. T(Q) C (or, equivalently, T(V + ) C V + and ( u o T)(u) ^ u(u>) for all u> £ V + ). We can consider the 
case where several transformations are applied because compositions of transformations are transformations again: If 
T, T' are linear maps V —> V which map fl inside then the same is true for the composition T o T' (we denote 
the composition of maps by a o symbol). 

We assume that the three parties act individually at spatially separated locations. Relativistic considerations lead to 
the consistency requirement that transformations performed by different parties must commute, e.g. if Alice performs 
a transformation Ta € Fa and Bob performs a transformation Tb £ Fb, then the total transformation must satisfy 
Ta°T b = T b o Ta- 

For our purposes, we do not need to specify the sets Fa, Fb and Fe any further; the only requirement is that 
transformations of distinct parties commute. The sets Fa, Fb and Fe define the systems A, B and E, i.e. we define 
the individual parties via the transformations that they can perform. This leads us to the following definition. 

Definition C . 6 : A tripartite scenario is a quadruplet 

Sabe = ((V 1 V + ,u),Fa,Fb,F e ), (C.9) 

where (V, V + ,u) is an abstract state space, and where 

Fa,Fb,Fe C {T : V -> V linear \ T(ft) C (C.10) 

are such that for all P, P' € {A, B , E} with P ^ P', it holds that Tp o Tp> = Tpi o Tp for all Tp £ Fp and for all 
Tpi £ Fp' ■ We call the elements of Fa, Fb and Fe the local transformations of A, B and E, respectively. 
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It is absolutely natural to define tripartite scenarios via commuting transformations rather than via a tensor product 
structure. In quantum theory, the two approaches are equivalent in finite dimensions (we will talk about this below). 
In more general infinite-dimensional cases, where it is not known whether the two approaches are equivalent, things 
are usually formalized in a commutative way rather than via tensor products (see [65], for example). Knowing about 
the equivalence in finite dimensions, we will formulate some quantum examples in the tensor product structure below. 

Example C.7 (A tripartite quantum scenario): One can formulate a tripartite situation in quantum theory as a 


tripartite scenario. Based on Example C.3, consider the tripartite scenario 

((Herm('H), Pos(?{), Tr), Ta, Tb,Te), where (C.ll) 

/ H='Ha®'Hb®'He , (C.12) 

Ta = {K-a <8> 1b <8> He | K.a is a trace non-increasing CPM on Herm('Hpi)} , (C.13) 

Tb — {1a <8>Kb <8> 1e I K b is a trace non-increasing CPM on Hermes)} , (C.14) 

Te = {1a ® 1b <8 K-e I Ke is a trace non-increasing CPM on Herm {'He)} , (C.15) 

where CPM stands for completely positive map. Having tensor product form, the local transformations of different 

parties commute. ■ 


For our purposes, Definition C.6 is all the structure one needs to specify. The local measurements are induced by 
the local measurements. We formalize this via the noation of a local instrument [56]. To get an intuition for what an 
instrument is, consider a Stern-Gerlach experiment. A spin-1/2 particle enters a magnet and undergoes one of two 
transformations: It either gets deflected upwards or downwards. Which of the two transformations it undergoes is 
determined probabilistically. Then it hits a screen, which reveals which of the two transformations the particle has 
undergone. This way, a measurement has been performed in two stages: a probabilistic application of a transformation 
and a detection. The sum of the probabilities of detecting the particle at the top or the bottom of the screen is one. 
If the state of the particle is described by a state oj G fi of an abstract state space, we may model this by a set of two 
transformations {P up , Tdown}- Such a set is an instrument. The norm u(T up (u>)) is the probability that the particle 
is deflected upwards, and likewise for u(T ( i OW n(w)). Thus, u can be seen to play the role of the screen, detecting the 
particle. The requirement that the particle must undergo one of the two deflections reads uoT up + it oXd own = u. The 
transformation T up is the analogue of the transformation p > P up pP up in quantum theory, where P up is the projector 
onto the spin-up state. Since u is given by the trace in quantum theory, the probability for the upward-deflection to 
occur is given by Tr(P up pP up ) = Tr(P up p), which is precisely the Born rule. 

A local instrument is such a set of transformations where all the transformations are the local transformations of 
one party. This motivates the following definition. 

Definition C.8 (Local instruments): For a tripartite scenario Sabe = {(V,V + ,u),T a,Tb,Te) with Q as defined 
in Definition C.2, we define the local instruments as the elements of 


Ip :=|lpCTp finite 


u o Tp = ■ 

Tp£lp 


Example C.9 (Local instruments in a tripartite quantum scenario): 
pie C. 7, we get that the local instruments are given by 


Dp = s Ip Q Tp 


T P is a TPCPM l 

TpCIp ) 


for P G {A,B,E} . (C.16) 

Considering the tripartite scenario of Exam- 

for P G {A, B,E}. m 


Remark C.10 (Local measurements): The definition of local instruments gives us a notion of local measurements 
as well. Consider a tripartite scenario Sabe = ((V, V + ,u),Ta,Tb,Te ) with its set of measurements M. It is easily 
verified that for a local transformation Ta &Ta, the map uoTa is an effect (as defined in Definition C.2). Likewise, for 
a local instrument I a G I a, the set {uoTa \ Ta £ I a} is a measurement. We interpret it as a measurement performed 
by Alice. We can also consider composite measurements where several parties locally perform measurements. For 
local instruments I a G Ia and Ib G Ib, for example, the set {uoTa° Tb \ Ta G I a, Tb G Ib} is a measurement. We 
interpret it as a composite measurement where Alice and Bob each perform local measurements, described by I a and 
Ib- The analogous holds for other parties and combinations thereof. 
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Example C.ll (Local measurements in a tripartite quantum scenario): Based on Examples C.7 and C.9, we can 
say how local measurements look like in a tripartite quantum scenario. A local effect of Alice is of the form 

Pabe | -t Tr(72.^4 ( g ) t B < S > Ie(pabe)) (C.17) 

for a trace non-increasing CPM IZa on Herm ('Ha)- However, for every such CPM, there is a POVM element Pa on 
Ha such that 6 


Tr {{Pa ® 1b < 8 > 1e)pabe) = Tx{JZa ( g > 1b < S > 1e{pabe )) • ( C . 18 ) 

This recovers the Born rule. Analogously, a composite measurement where Alice and Bob each perform local mea¬ 
surements consists of local effects of the form 

Pabe | -t Tr(Ki (g) TZ b <g Ie(pabe)) = Tr (( P ^ (g> Pg <8> 1e)pabe) ( C . 19 ) 

for POVM elements Pa, Pb on Ha, Hb- Thus, in our tripartite quantum example, local measurements reduce to 
POVM measurements of product form. ■ 

In Examples C . 7 , C.9 and C . ll , instead of choosing a tensor factorization for H and setting the local transformations 
to be acting non-trivially on one tensor factor, we could have chosen sets of transformations that merely commute, 
without a tensor product structure. The question of whether the resulting measurement statistics in that case would 
be different from the case with the tensor factor structure is known as Tsirelson’s problem [66, 67]. More precisely, 
the question is the following. Let H be a Hilbert space, let p be a density operator on H, let { Pk}k , {Qi}i be POVMs 
on H such that PkQi = QiPk for all k, l. Tsirelson’s problem is: Does there necessarily exist Hilbert spaces Ha, Hb, 
a density operator cr on Ha ®Hb and POVMs {Rk}k on Ha and {Si}i on Hb such that Tr (PkQip) = Tr{{Rk (g> Si)a) 
for all k, 17 In the case where H is finite-dimensional, the answer is known to be affirmative. For infinite-dimensional 
Hilbert spaces, the answer is still unknown. 

Thus, for finite-dimensional quantum systems, we can restrict ourselves to the case with the tensor product structure 
without loss of generality. For abstract state spaces, however, an analogous restriction might cause a loss of generality. 
The advantage of our weak definition of a tripartite scenario is that we do not need to know the answer to an equivalent 
of Tsirelson’s problem for generalized probabilistic theories. The downside is that it makes defining an equivalent of 
the min-entropy more difficult. We will deal with this issue in the next subsection. 

Notation: From now on, whenever we speak of a tripartite scenario Sabe, we implicitly assume that all its parts 
and induced structures are denoted as in Definitions C.l, C.2, C.6 and C.8 without restating it, i.e. instead of writing 
“Let Sabe = ((V, V + ,u),T a,T b,Te ) he a tripartite scenario, let fl be its set of normalized states, ... ”, we will only 
write “Let Sabe he a tripartite scenario”. 


2. A decoherence quantity for GPTs 

a. Motivation of an expression that quantifies decoherence 

We are now going to motivate an expression for the central quantitiy Dec(A\E) u for our decoherence analysis for 
GPTs. We take our inspiration from expression (A. 14) for the quantum min-entropy, which we repeat here for the 
reader’s convenience: 


H min (A\E) p = -\ogd A max F 2 {<h A A', 1a <S> 1Z e ^,A'{pae)) ■ (A. 14 revisited) 

F E _, A I 

There are two issues that prevent us from directly translating expression (A. 14) into our framework. The first issue 
is that in equation (C.8), to keep our framework as general as possible, we have defined a tripartite scenario with an 
overall state space (V, V + ,u) with tripartite states f2. We do not have notions of individual state spaces at hand. 
Thus, we do not have an analogue of a reduced state pae or of a transformation 7 Ze^A' from one state space to 
another. 


This can be seen from the Kraus representation of TZa- Tr (IZa(pa)) = Tr Ff,p a Cl) = Tr T 7 ) F,.p _\) = Tr (PaPa) f° r Pa = 
T., F^Fk- (We omitted the other tensor factors for brevity.) 
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The second issue is that we do not know what the analogue of a maximally entangled state $, 4 , 4 ' in our framework 
is. We resolve the first issue here in equation (C.19), arriving at an expression for Dec(A\E) w . In equation (C.22), we 
will then define what a maximally entangled state is in our framework. 

Expression (A. 14), which involves the state pae and TPCPMs IZe^A', can be transformed to an expression in 
which both the state and the TPCPMs are purified (see Figure 9). This expression will be our motivation for the 
expression for T)ec(A\E) u . The maximization over TPCPMs from E to A' is replaced by a maximization over unitaries 
from EE" to A' A", where E" and A" are ancilla systems extending system E and A', respectively. This is precisely 
the purification (or Stinespring dilation) of a channel as in equation (15). Since systems EE" and A'A" have the 
same dimension, we can identify their Hilbert spaces and regard the resulting Hilbert space as the Hilbert space of 
a system E t 0 t- This system involves all subsystems that the third party needs to control in order to bring itself as 
close as possible to maximal entanglement with Alice. Since UE tot is a transformation on system E tQ t alone, we can 
translate it into our generalized framework. 



A - >A A -> A 

(a) Unpurified situation (b) Purified situation 


FIG. 9: Purification of equation (A.14). Part (a) shows the system and maps involved in expression (A.14) for the quantum 
min-entropy. In expression (C.20), we purify this situation, as shown in (b), to arrive at a situation with three parties A, B 
and Etot, and with a map UE tot which acts on one system E to t alone rather than mapping from one system to another. 

The state pae is replaced by a purification pabe- We choose the purifying system B to be the channel’s output 
system, which gives us the overall picture of our decoherence analysis as shown in Figure 10. 



FIG. 10: Overall picture of our decoherence analysis for GPTs. Since the purifying system B in expression (C.20) is 
not specified, we can choose it such that it fits our situation for the decoherence analysis. 

The following gives a precise formulation of the purification of expression (A. 14). It can be proved using purification 
and Stinespring dilation. Let Ha, He be finite-dimensional Hilbert spaces of dimensions dA, dE, respectively, let 
Pae G S(Ha®He)- Then, for any purification pabe G S(Ha <8> Hb ® We) of pae, any Hilbert spaces Ha' , Id A" and 
He" of dimension dA' = dA, dA" = AaAe and dE" = dr A , respectively, any maximally entangled state <5>aa' G T aa' 
and any pure state |0)(0|s" G S(He"), it holds that 

H min (A\E) p = - log d A max max F 2 ($ A A' ® vba",(1ab ® U Etot )pABE tot (lAB ® Ut )) , (C.20) 

UE tot a BA" tot 

where pABE tot = Pabe ® |0)(0|s" and where the first maximization ranges over unitaries 

UE t ot '■ dd E ® He" —> H a' ® H a" , where Ha 1 ®Ha" —H E ®He" ='■ H Eto t (C. 21 ) 

and the second maximization ranges over pure states ctba" G S {Hb ® Ha" ) • 

Now we translate expression (C.20) into our generalized framework. We interpret the system E to t as the system 
controlled by Eve, and therefore rename E tot -A E. 

• Since we want to arrive at an expression that does not make unnecessary assumptions about the mathematical 
description of the physical situation, we avoid the factor dA present in (C.20). We look for a GPT analogue of 
max( 7 Etot max ffBJ „ E 2 ($ AA / ®<t B a", (Iab ® U Etot )pABE tot (Ias ® c4 tot )), omitting -logd^. As a consequence, 
we will have H m i n (A\E) p = — logdADec(A|£ , ) p in quantum theory (see Example C.14). 
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• We replace the maximization over all unitaries UE tot acting on system E tot by a supremum' over all local 
transformations Te £ Te- 7 8 

• We generalize the quantum fidelity to the fidelity in abstract state spaces as defined in Definition C.4. 

• We replace the state pABE tot = Pabe ® |0)(0|e" by a state u £ fi. 

• If we look at the state &aa' <£> &ba", we see that it is a state of maximal entanglement between Alice (A) and 

Eve (A'A") in the sense that by performing measurements with elements of the form Pa < 8 > Pa> ® Is ® 1 a" , they 
can get any statistics that two parties A and A' would be able to get by performing local measurements on the 
maximally entangled state <&aa'- We translate this into our framework by assuming that there is a set 4 'ae of 
“states with maximal correlation between Alice and Eve”. Instead of minimizing over states ® <Jba" , we 

then minimize over the set T ae- 

We postpone the discussion of how such a set ^ae looks like. We will give a definition of such a set in equation (C.22) 
below. For now, we write down an expression for our decoherence quantity Dec(A\E) u that depends on the choice of 
such a set C fl. According to what we have just discussed, the expression is 

sup sup F 2 (iP,Te(lo)) . (C.22) 

Te€l1~E 


We interpret the decoherence to be high when this quantity is high and vice versa, which is the opposite of H mm (A\E) p 
(see the end of equation (15)). Before we can define Dec(A|£%, however, we need to specify what a maximally 
entangled state in a GPT is. 


b. Definition of maximal correlation in GPTs 


The expression (C.22) for our decoherence quantity Dec(A|£’) w contains a maximization over a set ^ae C fl which 
we interpret to be the set of states with maximal correlation between Alice and Eve. We now define this set. 

Definition C.12: For a tripartite scenario Sabe, we define the set 'I 'ae of states with maximal correlation 
between Alice and Eve by 

For every binary local instrument I a = £ Ia, there is a binary local) ,,, . 

instrument Ie = {T%, T^} £Ie such that (uoT^oT^)(tp) + (uoT\oT^-)(ip) = 1 . J 

Definition C.12 can be read as follows. The superscripts 0 and 1 of the elements of the instruments I a and Ie stand 
for measurement outcomes, so (uoljo Tg)(ip) or (u o T\ o Tg)(ip) is the probability that Alice and Eve both get 
outcome 0 or both get outcome 1 , respectively, when they measure with respect to I a , Ie, respectively. Thus, the sum 
of these probabilities is the probability that Alice’s and Eve’s measurement outcomes are perfectly correlated. This 
means that for a state if £ ^ ae, it holds that for every binary measurement of Alice, there is a binary measurement 
for Eve such that their measurement outcomes are perfectly correlated. 

A closer look at some subtleties is advisable here, both to avoid confusion and to see the advantages of the weak 
assumptions that define our framework. With reference to Example C.ll, one may point out that that the set 

For every binary POVM { P {) A , P \} on Pa, there is a binary POVM 
{Pe, Pe} on He such that 

Tr((P 4 <g> 1 b <8> Pe)o) + Tr((P^ ® t B (g> Pe)o) = 1. 


£ S(Pa ® Pb ®Pe) 



'k AE := sip £. H 


7 We do not assume enough about e to guarantee that the maximum is achieved, so we replace it by a supremum. 

8 One might raise the objection that in the quantum case, Example C.7, the unitaries only correspond to those elements of 7~ e which 
bijectively map the space of density operators onto itself. It would be possible to include this restriction, but we decide not to do so, 
for two reason: We want to keep things simple, and we want to avoid the assumption that actions that the third party can perform can 
be purified as in the quantum case. 
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is empty. This may seem to make our definition of "F ae incompatible with quantum theory. Note, however, that the 
set 


| cr £ S(Ha <S> H b ® He) 


For every binary projective measurement {, P \} on Ha, there is a 
binary projective measurement {P^Pj,} on He such that 

Tr((P)} <8» 1b <8> Pe)<?) + Tr((P^ ®> ® Pj;)cr) = 1 . 


(C.25) 


is not empty as long as dim He ^ dim P 4 . If, as in equation (C.19), He = 'Ha' < 8 > Pa" with Pa' — Pa, then this set 
contains all the states of the form $aa' ® <tba" with <E>aa' G Faa' as in (C.20). The advantage of our weak definition 
of the local transformations is that it does not force to see Ta as the analogue of the set of all CPMs of the form 
TZa <8> Is ® 1 _e, but that it can be considered to be the analogue of all such CPMs which induce a functional of the 
form a 1 —»• Tr(Pcr), where P is a projector. Example C.ll can be modified accordingly (see Example C.14 below). 
This makes our definition of T 4 e compatible with quantum theory. 

With Definition C.12 at hand, we are finally ready to define the decoherence quantity. 


Definition C.13: Let Sabe be a tripartite scenario, let ui G Cl. We define the the decoherence quantity of uj by 

Dec(A|P) w := sup sup F 2 (ip, Te(lo)) (C.26) 

Te€l 1~E IpC.'&AE 


Example C.14: We consider a special case of a tripartite scenario in quantum theory. Consider 

((Herm(P), Pos(P),Tr), Ta,Tb,T e) , where 

H = Ha ® Hb <S> He 

TZa is a trace non-increasing CPM on Herm(PA) such that there is a) 
projector P A on Ha with Tr(P A pA) = T r(JZ A (PA)) for all p A € S(Ha) ) ’ 

and analogously for He a nd He- In addition, we assume for simplicity that Ha — He- In this case, 


Ta — \ Ha ® 1b ® He 


(C.27) 

(C.28) 

(C.29) 


'I'ae = \ <t G S(Ha < 8 > Hb < 8 > He) 


For every binary projective measurement {Pj(,P^} on Ha, there' 
is a binary projective measurement {Pjj, Pj,} on He such that 


= {$A 


E 1 


Tr((P°#l B ®Pg)< 7 )+Tr((P^ 
&B I $A£ G r AE , &B G S(Hb)} , 


l B ® Pe)<7 ) — 1 ■ 


(C.30) 


(C.31) 


where Tae is the set of maximally entangled states on S(Ha ® He) analogous to (A. 15). For a pure state pabe G 
S(Ha ® Hb ®He), this gives us 


Dec(A|P) p = max max max P 2 ($ab ® Pb, Ha® Hb® H-e(pabe)) 

IZe &AE a B 

= ma xF 2 (<b A E, Ha® He(pae)) 

Tce 

_. 1 r,-H min (A\E) 0 

d A 


(C.32) 

(C.33) 

(C.34) 


Hence, H min (A\E) p = -logeUDec(A|.E), : 


3. Bounds on the decoherence quantity for GPTs 

The goal of this subsection is to derive an upper bound on Dec(A\E) u in terms of the CHSH winning probability of 
Alice and Bob. This is a practically relevant bound: On the premise that the channel behaves identically in multiple 
uses and does not build up correlations between different uses (such a channel is said to be iid, for independent and 
identically distributed), this winning probability can be estimated through repeated measurements on Alice’s and 
Bob’s side. What we show is that this estimate in turn gives a bound on Dec(A|£%. In this section, we formulate 
this bound as a minimization problem which we solve and interpret in equation (C.89). 

In the following, we derive a lower bound on — logDec(A|P) w . We make the convention that — logO = 00 , where 00 
is a symbol for which we accept the inequality 00 ^ r for every real number r. This lower bound on — log Dec(A|P) w 
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then gives us an upper bound on Dec(A|i?) w . In a first step, we bound the fidelity-based quantity — logDec(yf|£ , ) tlJ 
by a trace distance-based quantity. This has the advantage that the resulting optimization problems which give us 
the bounds can be solved using linear programming. 

Proposition C.15: Let Sabe be a tripartite scenario, let ui & LI. Then 

— logDec(A|.E) w ^ inf inf D 2 (iP,Te(u>)) . (C.35) 

Te^lTe 


The following lemma is useful for the proof of Proposition C.15 below. 

Lemma C.16: For all x £ (0,1], it holds that — log(a; 2 ) ^ 2(1 — x ). 

Proof. We have that — log(;r 2 ) = —21og(a;), so the claim is equivalent to 

(x — 1) — log(x) > 0 for all x € (0,1]. (C.36) 

The functions F(x) = log(ai) and G(x) = x — 1 are differentiable on K>o- Thus, by the fundamental theorem of 

calculus, it holds that for all x £ M>o, 

F{x) = F(l) + j f(y)dy, G{x) = G{l) + J g(y)dy where f{y) = -^-F{y), g{y) = ^G{y), (C.37) 

so for all x £ ( 0 , 1 ], we have that 

(x - 1) - log(:r) = G{x) - F(x) = g(y) - f{y)dy = - 1 - d V > 0 • (C.38) 

<0 for all ( 0 , 1 ] 

This proves the claim. □ 

Proof of Proposition C.15. Since the right hand side of (C.35) is a finite real number, the inequality trivially holds if 
Dec(A|£’) tJ = 0 by the above convention. Thus, we assume in the following that Dec(A|B) tJ > 0. We have that 

— logDec(A|B) w = — log sup sup F 2 (i/j,Te(oj)) (C.39) 

Te£1~e i/'G^ae 

= -log( sup sup F(ip,T E (u}))) (C.40) 

\T e GTe tpe'fAE J 

For x £ (0,1], it holds that — logx 2 ^ 2(1 — x) (see Lemma C.16). Thus, since 

sup sup F(ip, T e {ui)) £ (0,1], (C.41) 

Te€1~e i/'G^ae 

we get that 


— logDec(A|£ , ) w ^ 2 1 — sup sup F(ip,T E (u>)) 

y t e oTe V’G'I'ae 

= inf jnf 2(1 - F(if>,T E (u))) 

Te€Te i/'G^ae 

= inf inf sup 2(1 — blip, T e (ui)\M)) . 
Te&Te V’G'I'ae mom 


(C.42) 

(C.43) 

(C.44) 


For the Bhattacharyya coefficient b and the total variation distance d, it has been shown [ 68 ] that for any two 
probability distributions distributions, it holds that 2(1 — 6 ) ^ d 2 . Since this is true in particular for the two 
probability distributions that the measurement M induces on the states if an d T e (u>), we get that 

— logDec(A|B) w ^ inf inf sup d 2 (if, Te(lo)\M) (C.45) 

TeoTe V'S'I'ae mom 

= inf inf D 2 (’iP,Te(uj)) , (C.46) 

Te€lTe i/’G^ae 

as claimed. □ 
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The idea that the fidelity and the trace distance are related is not new. In quantum theory, the Fuchs-van de 
Graaf inequalities (FvdG) relate the two quantities [69]. Inequality (C.35) is not completely analogous to the FvdG 
inequalities: It makes use of the logarithm in (C.35), which allows to apply classical relations that lead to a stronger 
bound than with the application of the FvdG inequalities. 

For the bounds that we are going to derive, the notion of a non-signalling distribution is central. Our bounds are 
essentially minimizations of functions over sets of non-signalling distributions Pr[a, b, c\x, y, z] u and Pr[a, c|x, z]^_, with 
certain additional properties. 


Definition C.17: A set of numbers Pr[a, b, c\x, y, z] u £ [0,1], indexed by numbers a,b,c £ {0,1} which we call 


outcomes, and numbers x,y,z £ {0, 1 } which we call settings, is a non-signalling distribution if 

normalization: E Pr t a , b, c\x, y, z] u = 1 for all x,y,z £ {0,1} , (C.47) 

a,b,c 

no-signalling: E Pr[a, 6 , c|0, y, z\uj y ^ Pr[(z, 6 , c| 1, y , zfor nil b , c, y, z G {0,1} , (0.48) 

a a 

^ Pr[a, 6 , c\x, 0, z] u = ^ Pr[a, 6 , c|x, 1, z]^ for all a, b,x,y G {0,1} , (C.49) 

b b 

E P r [ a , b, c\x, y, 0 ] w = E Pr t a ,b,c\x,y,l] u for all a, b, x, y £ {0,1} . (C.50) 

b b 

Similarly, a set of numbers Pr[a, c|x, z\^ £ [0,1], indexed by outcomes a,c £ {0,1} and settings x,z £ {0,1} is a 
non-signalling distribution if 

normalization: ^ Pr[a, c|x, z]^ = 1 for all x, z £ {0,1} , (C.51) 

a,c 

no-signalling: E Pr[a, c| 0 , z\$ = E P r [ a i c| 1 , z]$ for all c,z£ { 0 , 1 } , (C.52) 

a a 

EP r M*, 0 ]^=E Pr [a. c|a;, 1]^ for all a,x £ {0,1} . (C.53) 

C C 


The interpretation of equations (C.48) to (C.50) is that it is impossible for each of the three parties to signal to 
the other two parties by influencing their measurement statistics with the choice of the measurement setting. These 
one-party no-signalling constraints imply all the multi-party no-signalling constraints, saying that no collection of 
parties can signal to the remaining parties [70], so we do not need to require these constraints separately. 

Now we are going to formulate the bound on — log Dec(A|£% in terms of the CHSH winning probability of Alice 
and Bob. Assume that Alice, Bob and Eve are in a situation described by a tripartite scenario Sabe • Suppose that 
Alice and Bob have estimated that for the state ui £ fl that they are analyzing, their CHSH winning probability is at 
least A for some A £ [0,1]. Formulated in our tripartite scenario language, this means that they have found out that 
for local instruments 

I°A = {tT,Ta°} e Xa, I°b = {Tb 1 °, t b°} e X B , (C.54) 

I a = {tT’T 1 /} e Xa, lb = £ 1 B , (C.55) 

it holds that 

JE E (uoTfoT^){u)>\. (C.56) 

x,y a,b 
a(Bb=xy 

In that case, what can Alice and Bob infer about — logDec(A|A) w ? We have seen in Proposition C.15 that this 
quantity is lower bounded by ini TEe j- E inf^, e ^ 4B D 2 (iI>,Te{uj)). Alice’s and Bob’s estimate on their CHSH winning 
probability can be translated into a bound on this quantity. This is shown by the following proposition. 


Proposition C.18: Let Sabe be a tripartite scenario, let u £ fl be a state. If the CHSH winning probability of Alice 
and Bob is at least X, i.e. if there are local instruments /]}, I\, I B and I^ as in (C.54) and (C.55) and a X £ [0,1] 
such that (C.56) is satisfied, then 


inf inf D(iP,Te(oj)) ^ 
Te^'Te 


mm 

x,z£{ 0 , 1 } 
Y > Y[a,b,c\x,y,z] U j€:'D U) (A) 
Y>Y[a,c\x,z]ipeT>^ 




Pr[a, c\x, z - E Pr t a > c \ x > V> z b 

b 


(C.57) 
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where T> ul ( A) is the set of non-signalling distributions for Alice, Bob and Eve such that Alice and Bob have a CHSH 
winning probability of at least A, i.e. 




Pr[a, b , c|x, y, z]^ is a non-signalling distribution such that 

Z> W (A) = < 

Pv[a,b,c\x,y,z] u 

\22 E E Pr [ a ’ 6,c \ x ,y,z] u >\, 

x,y a,b c 



a(&b=xy 

> 


and where Vp, is the set of non-signalling distributions for Alice and Eve such that their measurement outcomes are 
perfectly correlated when they choose the same measurement setting, i.e. 


— 


|pr[a, c\x, z]^ 


Pr[a, c\x, z]^ is a non-signa.lling distribution 
such that Pr[a = c\x = z]^ = 1. 


(C.59) 


Proposition C.18 reduces our problem of lower bounding the decoherence quantity for GPTs to an optimization 
over non-signalling distributions. This allows us to use linear programming techniques, which in similar ways have 
been used in [26] to answer questions about non-signalling distributions. 

We need the following lemma for the proof of Proposition C.18 below. 


Lemma C.19: Let Sabe be a tripartite scenario, let upip £ f2. Then, for all local instruments ( Ia,Ib,Ie ) £ 
I A x Ib x Ie, it holds that 




Ta£Ia 

Ue^Ie 


(u O T a o U E )(tp) - 22 {uoT a °T b oU e o T e )(u) 
Tb£Ib 


(C.60) 


Proof. It is sufficient to show that for all u, ip £ fi, for all T E £ T E and for all (I a , Ib,Ie) £ I A x I E x I E , 


D{if,T E {u))> 


1 

2 


E 

Ta£Ia 

Ue^Ie 


(u o T a o U E ){xf) — 22 {uoT a oT b oU e oT e ){u) 
Tb£Ib 


(C.61) 


This is what we are going to show now. Let u,tp £ Cl, let T E £Te, let (Ia, Ib, Ie) G 1a x Ib x X e . Then 
D(ip,T E (u)) = sup d(ip,T E (w)\M) = )- sup 22 \e(ip) - e(T E (w))\ . 


MeM 


MeM 


(C.62) 


eSM 


If instead of taking the supremum over A4, we only evaluate the expression for a particular element of A4, we get a 
lower bound on (C.62). We choose the element (c.f. Remark C.10) 


Hence, 


{« o 1 a ° U E | Ta £ Ta, U E £ 1 E } £ M . 


D(ip, T e (u)) ^ - 22 \{u°T A ° U E )(tp) - {uoT a o U e oT e )(u)\ . 

Ta €Xa 

u E ex E 


By the definition of a local instrument, u = )T) r eI u o T b . Thus, 


DW,T e (w))>- 22 


Ta €Xa 
U e^-Te 


(u o T a o U E )(yp) — 22 ( u ° T B o T a o U e o T e )(u) 
Tb€Ib 


(u O Ta o U e )(4>) - 22 ( u °Ta°T b oU e oT e )(u) 

Tb£Ib 


-I £ 

Ta (zT/ 

u E eXi 

where in the last equality, we made use of the fact that transformations of different parties commute. 


(C.63) 


(C.64) 


(C.65) 

(C.66) 

□ 
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Proof of Proposition C.18. It is sufficient to show that for every xf> £ 5 'ae, the claimed inequality holds without the 
minimization over ^ ae, he. 


inf D(tjj,T E (u ;)) ^ 

Te€lT e 


mm 

x,zE{ 0,1} 

Pr[a,b,c\x,y,z] UJ €:'D U j (\) 
Pr[a y c\x,z] lJj eT>^(^AEA) 


It 


Pr[a 


Z U — Pi'[a,, 6, 


c\x,y,z\ u 


(C.67) 


By means of Lemma C.19, we know that for all x, y £ {0,1} and every I E £ P E , 


1 


iffi^W)> 2T mf„ E J2 


Te^lT e 


aG{0,l} 

Ue^Ie 


(«olfo(/ £ )(^)- Y, (uoTf oT b J y oU e oT e )(co) 

be{ 0,1} 


Let V’ € ^ae, let I E = {Ug°, I E = {Ug 1 , U^ 1 } £ I E be local instruments for Eve such that 


{u o tT o u°j°m + {u o rf o ufm = i, 

(u O Tf o ufm + (u O Tf o ufm = 1, 


L A ^E 

rfoUf)^) + {uoTfoUf 

which exist according to the definition of ^ae (Definition C.12). It holds that for every x,y,z £ {0,1}, 


^ \teGTe E 


Te^lT e 


1 


a,cG{0,l} 


= - inf V 
2 T E er E ^ 


a,cG{0,l} 


(uoT% x oU% z ){if)- Y (« oT i oT b oUf oT e )[lo) 
&e{o,i} 


Pr[a, c\x, z\$ - Y Pr [a,b,c\x,y,zl 
&6{0,1} 


where 


Pr[a, c\x, z]^ = (mo ^ 1 o U$ z )(ip ), 
Pr[a,6,c|a:,2/,z] w = {u o T^ x o T^ v o U^ z o T E )(u). 


(C.68) 


(C.69) 

(C.70) 


(C.71) 

(C.72) 


(C.73) 

(C.74) 


Hence, 


inf 

TsGT e 


D{ip, T e (u))) 




— min inf 

2 x,y,z T E eT e 


E 

a,cG{0,l} 


Pr[a, c\x, z\^ 


Y ^[a,b,c\x,y,z) u 
6e{o,i} 


(C.75) 


Pr[a, c\x, z\f forms a non-signalling distribution: For the normalization, note that for all x, z £ {0,1}, we have that 


E Pr t a ’ c \ x > = E( u ° t a x ° u e z )w 

a,c a,c 



Y uoU e z (VO = W (V0 = 1 • 


For the no-signalling condition, note that for all c,z£ {0,1}, it holds that 


Y] Pr[a, c|0, z\^ 


((E uoT a | 0 )°^ wo 

^Pr[a,c|i,2]^ , 


( uoU E Z )W 



(C.76) 


(C.77) 


(C.78) 


(C.79) 
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and that for all a,x £ { 0 , 1 }, it holds that 


E Pr t a ’ c \ x ' °iv> = E( u ° t a x ° u e°)w = E ( m ° u e°° t a x )w = ^|e u ° u E°^j ° t a x ^J wo 
= (uot“ |x )(v>) = ^E^ 0 * 7 ^ ° T A |x j (vo = 5 Z Pr [ a ’ c i ;Ei1 ^’ 


(C.80) 

(C.81) 


where in the second equality, we made use of the fact that local transformations of different parties commute. Anal¬ 
ogously, for every Te £Te, Pr[a, b, c\x, y, z] u is a non-signalling distribution. Moreover, Pr[a, b, c\x, y, z]^ satisfies 


^E E E Pr [ a ’ 6 ’ c l x ’V’ z \“ = \Y1 E '52( uoT a X oT b V oU e Z oT e){u) (C.82) 

x,y a,b c x,y a,b c 

a(&b=xy a(&b=xy 

= lE E IJ 2 ^ oU e Zo T E ) 0 Ta X °Tb V ^) ( C - 83 ) 

x,y a,b \ c / 

a^b—xy s * "v- y 

u 

3s A, (C.84) 

where the inequality is one of the assumptions of the proposition. Furthermore, Pr[a, c\x, 1]^ satisfies 

Pr[0,0|0,0 ] 0 + Pr[l, 1|0,0]* = (u o T°'° o V$ 0 )(if>) + (u o T 1 } 0 o = 1, (C.85) 

Pr[0,0| 1, 1U + Pr[l, 1|1, IU = (« ° Tf o T 0 ' 1 )^) + (u o T 1 / o uf)(if) = 1 (C.86) 

(where we made use of (C.69) and (C.70)), which we may abbreviate as Pr[a = c\x = z]^ = 1. Thus, for every 

Te £Te,w£ have that Pr[a, b, c\x, y, z\ u € V bJ ( A) and Pr[a, c\x, z]^ € Thus, 


inf 

TeGT e 


D(t/>,T E (u)) 


> 


— min inf 

2 x,y,z Pr[a,6,c|a:,2/,z]„eDu(A) 
Pr[a,c\x,z] 1 p£'D,j,('i’ ae,\) 


E 

a,cG{ 0 ,l} 


Pr[a, c\x, z\ip — E Pr[a,b,c\x,y,z] u 
be{o,i} 


(C.87) 


Since Pr[a, b, c\x, y, z] u satisfies the no-signalling property, the right hand side of (C.87) is independent of y, so the 
minimization only needs to be performed over x and z. Moreover, the infimum over the sets T> u ( A) and T>^ is a 
minimum because it is the infimum of a continuous function over a convex polytope, which is always attained (see 
equation (C.89) for more details). This completes the proof. □ 

Corollary C.20 (The bound): Let Sabe be a tripartite scenario, let oj £ fl be a state. If the CHSH winning 
probability of Alice and Bob is at least A (in the above sense), then 

Dec{A\E) u < 2~ s2 ^ , (C. 88 ) 


where 


(5(A) 


min 

x,z£{ 0 , 1 } 

Pr[a,6,c|x,y,2:] w G'D u; (A) 
Pr[a,c|:r,;z] 1 /,G'D</, 



a,c 


Pr[a, c\x, z}^ - J2^b,c\x,y,z\J , 

b 


(C.89) 


Proof. This is a direct consequence of Propositions C.15 and C.18. 


□ 


4. Evaluation of the bound and results 


a. Formulation of the bound as a linear program 

In this subsection, we evaluate the bound (C. 88 ). To this end, we rewrite (C.89) in terms of linear programs. 
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Linear Program: The bound 6(X), which is a function 5 : [0,1] —► [0,1], is given as follows. For all A G [0,1], the 
value 5(A) is the solution of the linear program 


minimize 
subject to 


S( A) 

Pr[o, b , c\x, y, z] u G T> u ( A) 

Pr[a,cja;, 2 ]^ G 

> Ea.c^ac Vz,z G {0,1} 

s ac > \ ( Pr K c \ x -> Z U - E b Pl 'k c \ x i 0, z\ u ) ^ ~5 XZ Va, c, + 2 : G {0,1} 


This is a linear program in 97 variables: 

{ Pr[o, b, c\x,y, z ]ui} abcxyze{01} 
{Pr[o,c|a:,|/]^} OiC>XiZe{0jl} 

{Ec } a,c,x,z£{ 0,1} 

5(A) 


64 variables 
+ 16 variables 
+ 16 variables 

+ 1 variable 
= 97 variables 


(C.90) 


We have already written out the constraints for these 97 variables as (in) equalities. The third and fourth line are 
already written as such in the program description, and for the first two lines, we refer to the following: 

constraint (in) equalities 

Pr[o, 6 , c\x, y, z] u G V U (X) (C.47) to (C.50) and inequality in (C.58) 

Pr[a, c\x, z\^ G V^{\) (C.51) to (C.53) and equation in (C.59) 

The inequalities define a convex polytope over which the convex function 5(A) is minimized, so the minimum is 
attained. It is straightforward to bring these inequalities into the standard form of linear programming. We solved 
the resulting linear program using standard linear programming routines in Mathematica and Octave. 


b. Solution of the linear program and discussion of the results 



FIG. 11: Plot of the result. The bound 2 5 ^ is non-trivial precisely when the CHSH winning probability for Alice and 
Bob is non-classical, i.e. A > 3/4. 


We plot the result in Figure 11. The bound 2~ s ^ is non-trivial for values A G (3/4,1]. This is a very satisfactory 
result as one cannot expect the bound to be non-trivial for A G [0,3/4]: A CHSH winning probability of at least 
A G [0,3/4] for Alice and Bob is always compatible with Dec(A|£') w = 1. To see this, note that the requirement 
for a state to yield a CHSH winning probability for Alice and Bob of at least A G [0,3/4] is trivial: Alice and Bob 
can choose trivial measurements that always yields 1 as an outcome, independently of the state. More precisely, in 
our tripartite scenario language, we can express this as follows. Certainly, there are tripartite scenarios in which the 
identity map ly and the zero map 0y are in Ta, Tb and Te - 9 For such tripartite scenarios, the condition (c.f. 


9 Every tripartite scenario can be turned into such by adding 1 y and Oy to the sets of local transformations. In fact, it would be physical 
to assume that each set of local transformations contains 11 y and 0y. 
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(C.54) to (C.56) in equation (C.34)) 


x,y 


,o|o 


nO|l rpllll 


E E ( u ° T a X o T b V ') [O’ 3/4] for some (C.91) 

a,b 
b—xy 


{t7,Tb 




is always satisfied because for all w G fl 
1 


j £ £ 

x,y a,b 
a(Bb=xy 


. . •» ri rO 10 /t-t0|1 

fo * T i\o~ A 


^ zj-iUjU _ _ 2~i0|0 _ rj~iyj\.L _ q 


- 0|1 _ ( 
B 


T'- l l u _ Til|l _ rp 110 _ rpl\l _ ^ 

1 A ~ 1 A ~ 1 B ~ 1 B ~ U-V ■ 


(C.92) 


This means that the requirement that the CHSH winning probability for Alice and Bob is at least A G [0,3/4] does 
not exclude the case u> € 'h ae • In that case, Dec(A|S) w = sup Te67 - e sup, 0 g ^, AE F 2 {iP,Te{cj)) = 1. 


Appendix D: A test for gravitational decoherence 
1. An optomechanical setting and its model for gravitational decoherence 

The objective here is to create two entangled photonic qubits in which one photon is prepared in an opto-mechanical 
system that is itself subject to gravitational decoherence — if there is any — and the other photon is prepared in 
an identical cavity except the mirrors are fixed and cannot move. This model is a modification of the model first 
proposed by Bouwmeester [5] in which an itinerant single photon pulse is injected into a cavity rather than created 
intra-cavity as here. Our modification avoids the problem that the time over which the photons interact with the 
mechanical element is stochastic and determined by the random times at which the photons enter and exit the cavity 
through an end mirror. In the new scheme, the cavities are assumed to have almost perfect mirrors — very narrow 
line width (see for example [35]). 

The intracavity single photon Raman source is described in Nisbet-Jones, et al. (34]. In this scheme (see Fig. 12) a 
control pulse can quickly and efficiently prepare a cavity mode in a single photon state by driving a Raman transition 
between two hyperfine levels we label as |g), |e). In our scheme there are two optical cavities otherwise identical except 
in one of the cavities a mechanical element can respond to the radiation pressure force of light. 

We will assume that we can prepare the atomic sources in an arbitrary entangled state \g, e) + | e,g), for example, 
using the trapped ion schemes of Monroe [71]. In addition we will assume that we can make arbitrary rotations in the 
g , e subspace of each source and also make fast efficient single shot readout of the state of each source, for example 
using fluorescence shelving. This means we can readout the atomic qubit in each cavity in any basis. 

The write laser implements the Hamiltonian H w = *MI(f)(aI|e)(c/| — a\g)(e\)/2. This is a rotation in the state space 
{\g)\0), |e)11)}. We can thus prepare arbitrary states of the form cos0/2|5 , )|O) + sin0/2|e)|l), where 9 is determined 
by the pulse area. We will refer to the case of 9 = n as a 7r-pulse. Note that if the source is in the excited state |e) 
and the cavity is in the vacuum, no photon is excited. 

Starting with the cavities in the vacuum state the protocol proceeds as follows: 

1 . prepare the source atoms in the state | g,e) + \e,g). 

2. apply the write laser with a 7r-pulse 

3. free evolution of the OM systems for a time T 

4. apply the write laser with a 7r-pulse 

5. readout the atomic state in each cavity. 

At the end of Step 2, the state of the sources and the cavities is ^ 2 ) = |e, e) (g> (| 1,0) + |0,1)) where |n, in) = \n) (g> |m) 
with each factor being a photon number eigenstate. 

a. Gravitational decoherence. 

We will use Diosi’s theory of gravitational decoherence[l]. This is equivalent to the decoherence model introduced 
in Kafri et al. [9]. One mirror of the opto-mechanical cavity is free to move in a harmonic potential with frequency 
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Raman single photon source 



E(t) 


a t 1 


-9 / 


A- 1 ~ 1 

—-*- • 



FIG. 12: Two cavities each contain a Raman single photon source controlled by an external laser ‘write field’ E(t). The Raman 
sources are first prepared in an entangled state. Only one cavity contains a mechanical element coupled by radiation rouser to 
the cavity field. 


u! m . The master equation for a massive particle moving in a harmonic potential, including gravitational decoherence 
is 


— = p\ - A grav [6 + b\ [6 + b\p]\ 


(D.l) 


where 


b = 


ITlIjJm 

2 h 


x + i 


\/2h\ 


mujn 


=P 


(D.2) 


with x, p the usual canonical position and momentum operators. The gravitational decoherence rate A grav is given 

by 


, 2ttGA ^ N 

Ag ra v = ~X~ - (D.3) 

O id-m 

with G the Newton gravitational constant and A the density of the mechanical element. As one might expect A grav 
is quite small, of the order of 10 -8 s^ 1 for suspended mirrors ( as in LIGO) with w m ~ 1. 

Form a phenomenological perspective the effect of gravitational decoherence is analogous to a Browning heating 
effect. To see this we note that the average vibrational quantum number increases diffusively 

= 2A grav (D.4) 

Indeed, one could simulate this effect by adding a stochastic driving force to the mechanical element via the stochastic 
Hamiltonian 

Hs = f (6 + 6t) (D - 5) 

where I{t) satisfies an Ito stochastic differential equation, 

dl(t) = \/4A^A dW(t) (D.6) 

where dW(t) is the Weiner increment. Averaging over all histories of the stochastic driving force gives the final term 
in Eq. D.l. 
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In the absence of mechanical dissipation, there is no steady state. In reality the mechanical quality factor, Q = 
Wm/lm, is finite leading to a steady state with mean phonon number given by 

<& f b) ss = (D.7) 

'Jm 

This of course assumes that there is no additional mechanical heating (regular thermodynamic kind): hardly a realistic 
assumption. This adds a large (comparatively) additional term to A grav so that we find (for >> huj m ), 

ksT 

Agrav t Ag rav T Aheat i where A heat, ~ }l(f * (D.8) 

Given the incredibly large quality factor of Q = 10 10 , one would need to cool the mechanical element to nano-Kelvin 
for the thermodynamical heating to be of the order of the gravitational heating. 


b. Optomechanical probe of gravitational decoherence. 

The optomechanical Hamiltonian in cavity-one is 

Hom = hjj m tfb + hgo{b + tf) (D.9) 

go is the single photon optomechanical coupling rate. Typically go ~ 1 s -1 for the sorts of cavities we are considering 
here. This is about the same order of magnitude as ui rn . In new field OM cavity technologies, go can be as high as 
10 3 s _1 however in such cases the mechanical frequency is also typically much higher ~ tens of MHz. The interaction 
time is T which is short compared to the cavity decay time (which we neglect). We will assume that the mechanics 
starts in a thermal state, the steady state of the system subject to gravitational decoherence, heating and dissipation. 
This is given by 


1 f . a | 2 

Pom = — d 2 a e~~*~ |a)b(a| (D.10) 

t rn J 

where n = (b^b) ss is the steady state mean phonon number given in Eq. D.7. 

It is simplest to work in an interaction picture defined by the mechanical free dynamics, 

Hom,i = hgoibe-^ + (D.ll) 

The corresponding unitary evolution operator is 

U(t) = e mbi -P* (t)b (D.12) 

where 

/3(t) = — 1) (D.13) 

The initial state for the OM interaction is the state at the end of Step 2 

Pom(0) = x - (|1, 0)(1, 0| + |0, 1)(0, 1| + |1, 0)(0, 1| + |0, 1)(1, 0|) ® Pm (D.14) 

where p. m is the state of the mechanical element at the start of the protocol, a thermal state. We can ignore the state 
of the atomic sources at this stage as they do not participate in the OM interaction. 

The state of the optomechanical system after an interaction time T is given by 

Pom(t) = i(|l,0)(l,0|p m + |0,l)(0,l|t/(t)p m 17 t (t) + |l,0)(0,l|p m H t (t) + |0,l)(l,0|t7(t)p m ) (D.15) 

The reduced state of the cavity fields is given by tracing out the mechanical degree of freedom, 

Pf(t) = l - (|1,0)(1,0| + |0,1)(0,1| + i?*|l, 0)(0,1| + -R|0,1)(1,0|) 


(D.16) 
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where 

R = e -(l+2fi)|/3(t)| 2 /2 


where 


with (as above) 
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Continuing with the protocol from Step 4, now results in the state of the atom-field system 


Paf{t ) = Pa(t) <S> |00)(00| 


(D.17) 

(D.18) 

(D.19) 

(DAO) 

(D.21) 

(D.22) 


where 


Pa(t) = ^(\ge)(ge\ + \eg)(eg\ + R*\eg)(ge\ +R\ge)(eg\) (D.23) 

The suppression of coherence due to the thermal state of the mechanics has been transferred to a reduction of 
entanglement in the atomic sources. A readout of the atomic sources will reveal this through either state tomography 
or via a reduction in a CHSH correlation for a Bell-type experiment. 

The function R(t) is a periodic function of time. At each period of the motion it returns to its initial value of zero 
and the cavity field state would return to the fully entangled state it was in after Step 2. If we chose T = Zit/uim 
then the protocol will return the atomic system to the same entangled state in which it began. This is because we 
have ignored the heating of the mechanics over the period T so the only way decoherence enters is through the initial 
thermal excitation of the mechanics. In effect the protocol is a thermometer. We thus see that for maximum effect 
we need to ensure go u> m - On the other hand, gravitational heating requires a small value of uj m and typically 
such OM systems have go/uj m <C 1. Perhaps technical advances will enable OM systems with long mechanical periods 
and large single photon coupling. Of course this will also require sub hertz cavity line widths. In the (exceptionally) 
optimistic case we can take T ~ 1 nK, u m ~ 1 s -1 , Am ~ 10 -10 s -1 so that Q ~ 10 10 . 


2. An experimental test of the model 

In equation (C.92) above, an optomechanical setting has been described. Making some assumptions about how 
gravitational decoherence influences the optomechanical system, a model has been given that describes how the 
state of the optomechanical system changes over time. In this section, we consider this model for the state of the 
optomechanical system as given and analyze it using our decoherence test formalism. We calculate the amount of 
decoherence that would be introduced to the optomechanical system if the model was correct. We compare this 
to the amount of decoherence that one would observe if there was no such gravitational decoherence, determining 
the difference between the two predictions. We devise an experiment that aims at estimating the actual amount of 
decoherence at a point in time when this difference is maximal. This turns the optomechanical experiment into a test 
that allows to falsify the above model for gravitational decoherence if it was wrong. This shows that the decoherence 
testing formalsim presented in this work can be applied in situations where the physical process is unknown. It allows 
to subject proposed models of the process to a consistency check. 

We first present the predicted values of T)ec(A\E) p of the optomechanical system for the two cases where grav¬ 
itational decoherence is present or absent, respectively, for some example parameters of the experiment. We then 
calculate the CHSH value /3 that one would have to measure in order to falsify the model for gravitational decoherence. 
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The main quantity of interest in our analysis is the decoherence quantity Dec(A|£?) p for the state pab = p/(f) 
desribed in equation (D.16). This can be calculated using Lemma B.9. It turns out to be 

Dec{A\E) p = i (l + y/l - R^j (D.24) 

= \ ^1 + ^1 - exp ^—4(1 + 2n)^§- sin 2 • (D.25) 

If the above model is correct and gravitational decoherence occurs, both the gravitational interaction and the me¬ 
chanical heating contribute to the average vibrational quantum number, i.e. we have 


^ — r/.g ra v + Uheat 


47 tG 1 . 2k B l 

-w- -A+-A — 

*3 7m^m 


(D.26) 

(D.27) 


and thus 


Dec(A|£?) p = - 1 



1 — exp —41 + 2 



V 3 7mW. 


(gr. dec. included). 

(D.28) 


If gravitational decoherence is absent, then only the mechanical heating contributes to the average vibrational quantum 
number, i.e. we have 


71 — ^heat 


2 k B 1 

h U! m 


T 


(D.29) 

(D.30) 


and thus 


Dec(A\E) p = i ^1 + y 1 - exp ^-4 ^1 + 2 ~ T ^j ) ~ 2 ~ sin2 ) j (S r ‘ dec - neglected). (D.31) 

Figure 13 shows how the decoherence quantity in equation (D.28) as a function of time varies for different materials of 
the mechanical element and different temperatures, compared to the case where there is no gravitational decoherence 
as in equation (D.31). 

In order to rule out the model for gravitational decoherence, one needs to measure a CHSH value /3 which is 
incompatible with the value of Dec(A\E) p given in (D.28). The minimal value /?f a i s of /3 that needs to be measured 
for this falsification can be calculated using Theorem B.l: Using MATLAB, we numerically evaluated the quantum 
bound on Dec(A|£l) p , which is given as a point-wise maximization problem in Theorem B.l. We inverted the resulting 
set of data points and interpolated a function from the resulting data using Mathematica. The resulting function 
takes a value of Dec(A|i?) p as its input and outputs the minimal /? that needs to be exceeded in a measurement in 
order to rule out the given value of Dec(A|i?) p . Thus, applying this function to the curves of the Dec(A|i?) p values 
of the gravitational decoherence model in Figure 13 yields the curves for (3f a \s- The results are plotted in Figure 14 
for the same materials and temperatures as above. 

In order to determine whether it is promising to measure a value of /? that lies above /3f a i s , we need to determine 
the value /3 mec h of (3 which is predicted in the case where gravitational decoherence is absent, i.e. where we only have 
mechanical heating. We can do that exactly: Equation (D.16) gives us an expression for the state, which we consider 
for the value of R given in the case of mechanical heating only, n = n m ech- Then we calculate the value of /3 mec h for 
the case where the measurements are taken to be the standard CHSH measurements 


Aq — (JX , 


B 0 = 


&x — O’* 


Ai = a z , 


B 1 


&x + 

~^fT~ 


where cr x , a z are the Pauli x- and z-operator, respectively. The resulting curves are shown in Figure 14 as solid curves. 
It turns out that for the relevant time interval (where /3 mec h is larger than either of the /3f a i s ), the curve of /3 mec h for 
the standard CHSH measurements is almost identical to the curve one would get for the optimal measurements for 
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T = InK T = 5nK T = 10 nK T =< 50 nK 






-Gravitational decoherence included, density A = 2,102 x 10 4 kgm ~' 5 (rhenium) 

-Gravitational decoherence included, density A = 2, 7 x 10 s kgm -4 (aluminum) 

-Gravitational decoherence neglected, (material-independent mechanical heating only) 


FIG. 13: Predicted values of the decoherence quantity in the optomechanical experiment. The decoherence quantity 
r>ec(A\E) p as in Equation (D.28) is plotted as a function of time for different temperatures and two different materials of the 
mechanical element. In addition, Dec(A|j5) p is plotted for the case where there is no gravitational decoherence, Equation 
(D.31). The calculations have been made for the experimental parameters g 0 = 1 s~ 4 , cn m = 1 s _1 and 7 m = 10 -10 s _1 . 


each time t. The latter can be calculated using a formula presented in [24]. This is an experimentally desirable fact: 
Using a fixed measurement independent of the measurement time is almost optimal. 

The most promising measurement time for a falsification of the gravitational decoherence model is given by the 
time when /3 m ech (that one may hope to actually measure) is high but /3f a i s (which one needs to exceed) is low. Thus, 
the optimal measurement time can be calculated as the time t max that maximizes the gap function 

g{t) := /3 mec h(f) - /Jfais (t) ■ (D.32) 

This gap function is depends on the density A of the mechanical element and its temperature T. One can see that 
temperatures that look promising for a falsification measurement when looking at the T)ec(A\E) p values in Figure 13 
turn out to be too warm when looking at the experimentally relevant analysis of the /? values in Figure 14. As an 
example, we have calculated the optimal measurement times for T = 1 nK for the densities of aluminum and rhenium. 
They are visualized in Figure 15. If there is no gravitational decoherence, one needs to measure values of /3 that are 
~ 0.1 close (aluminum) or ~ 0.2 close (rhenium) to the value that one can maximally measure using the standard 
CHSH measurements, in order to exclude gravitational decoherence. 
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T = 1 nK T = 5 nK T = 10 nK T = 50 nK 






heat 0 ) 

- /3f a i s (t), density A = 2, 7 x 10 3 kgm -3 (aluminum) 

- /3f a is(t), density A = 2,102 x 10 4 kgm~ 3 (rhenium) 


FIG. 14: Minimal CHSH values for the falsification of the gravitational decoherence model. The quantity /3f a i a , 
which is the minimal value that needs to be exceeded in the measurement of the CHSH value /3 in order to rule out the 
gravitational decoherence model, is plotted as a function of time for the same materials and temperatures as above. In addition, 
the value /3 mec h is plotted, which is the CHSH value that can actually be measured using the standard CHSH measurement in 
the case where gravitational decoherence is absent and only mechanical heating contributes to the decoherence. 


T=lnK 



t 

FIG. 15: Optimal measurement times for ruling out the gravitational decoherence model. The three plots are 
identical to the ones in the leftmost box in Figure 14, i.e. for T — 1 nK. In addition, the time f max at which the gap g(t) 
between /5 me ch and /Jf a i a is maximal is indicated for the two cases where the material of the mechanical element has the density 
of aluminum or rhenium. 



















